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Abstract 

In the so called lightbulb process, on days r = 1, . . . , n, out of n lightbulbs, all initially off, exactly r 
bulbs, selected uniformly and independent of the past, have their status changed from off to on, or vice 
versa. With X the number of bulbs on at the terminal time n, an even integer, and fi = n/2, a 2 = Var(X) , 
we have 
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where Z is a standard normal random variable, and 
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yielding a bound of order 0(n _1//2 ) as n — > oo. A similar, though slightly larger bound holds for n odd. 
The results are shown using a version of Stein's method for bounded, monotone size bias couplings. The 
argument for even n depends on the construction of a variable X s on the same space as X that has the 
X-size bias distribution, that is, that satisfies 

E[Xg(X)] = [iE[g(X s )] for all bounded continuous g, 

and for which there exists a B > 0, in this case 5 = 2, such that X < X s < X + B almost surely. The 
argument for n odd is similar to that for n even, but one first couples X closely to V , a symmetrized 
version of X, for which a size bias coupling of V to V s can proceed as in the even case. In both the 
even and odd cases, the crucial calculation of the variance of a conditional expectation requires detailed 
information on the spectral decomposition of the lightbulb chain. 

1 Introduction 

The problem we consider here arises from a study in the pharmaceutical industry on the effects of dermal 
patches designed to activate targeted receptors. An active receptor will become inactive, and an inactive 
one active, if it receives a dose of medicine released from the dermal patch. Let the number of receptors, all 
initially inactive, be denoted by n. On each day of the study, some number of randomly selected receptors 
will each receive one dose of medicine, changing their statuses between the inactive and active states. We 
adopt the following, somewhat more colorful, though equivalent, 'lightbulb process' formulation from [7]. 
Consider n toggle switches, each connected to a lightbulb, all of which are initially off. Pressing the toggle 
switch connected to a bulb changes its status from off to on and vice versa. The problem of determining the 
properties of X, the number of light bulbs on at the end of day n, was first considered in [7J for the case 
where on each day r = 1, . . . , n, exactly r of the n switches are randomly pressed. 

More generally, consider the lightbulb process on n bulbs with some number k of stages, where s r € 
{0, . . . , n} lightbulbs are toggled in stage r, for r = 1, . . . , k; we refer to the vector s = (si, . . . , Sk) recording 



°AMS 2010 subject classifications: Primary: 62E17, 60C05 Secondary: 60B15, 62P10. 
°Key words and phrases: Normal approximation, toggle switch, Stein's method, size biasing. 



1 



the number of bulbs affected on each study day as the 'switch pattern.' In order to consider quantities that 
depend on some subset of size b of the n bulbs, we define 



X n ,b,s — 



t=0 



t) { ' (»), 



and A„ :bjS 



Y\ ^n,b,s r , 
r=l 



(1) 



where (n)k = n(n— 1) • ■ ■ (n — fc + 1) denotes the falling factorial, and the empty product is 1. Generalizing the 
results in [7] , writing X s for the number of bulbs on at the terminal time when applying the switch pattern 
s = (si, . . . , s n ), the martingale method in Proposition 4 of [TT] shows that if the process is initialized with 
all bulbs off, then 



EX S — — (1 — ^n,l,s) 



and Var(X s ) = J(l - A„, 2 , s ) + j(A n , 2 , s - A£ M ), 



where, from (fl}, 



2s 4s 4s(s - I) 

— and X n .2, s = 1 1 t tt- 
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for s = I. 
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Letting n = (1, . . . , n), we call the standard lightbulb process the one where s = n, and in this case 
we will write X short for X n . In particular, @ with s = n recovers the mean /i = EX and variance of 
a 1 = Var(X) as computed in [7J. Other results in [7] include recursions for determining the exact finite 
sample distribution of X. Though computational approximations to the distribution of X, including by the 
normal, were also considered in [7J, the quality of such approximations, and the asymptotic normality of X 
was left open. 

Theorem 1 1 . 1 1 settles the matter of the asymptotic distribution of X by providing a bound to the normal 
which holds for all finite n, and which tends to zero at the rate n -1 / 2 as n tends to infinity. We consider 
the case of even and odd n separately. In the even case we directly couple the variable X to one having the 
X-size bias distribution, as described later on in this section. In the even case © yields A Illn /2 = 0, and 
therefore A ni i jn = 0, hence from @ we obtain EX = n/2, and also that a 1 = Var(X) is given by 



n n 

— (1 — A„^2,nJ + — A n ,2,n- 



(4) 



To state our result for n odd, for s = (si, . . . , Sk) let 



. . r 



i) re {m, m + 1} 
X n ,b,r otherwise 
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that is, A nj f,. s is obtained from X n ,b,s by replacing any occurrences of X n ,b.m and X n- b, m +i in the product (TT]) 
by their average. In the odd case, we proceed by first coupling X to a more symmetric random variable V 
with mean and variance given respectively by 



fj, v = - and 



= 2(1 

4 V 



X, 
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Then, with V in hand, we couple V to a variable with the y-size bias distribution, and proceed as in the 
even case. In Theorem ll.il and the remainder of the paper, Z denotes a standard normal random variable. 
Theorem 1.1 Let X be the number of bulbs on at the terminal time n in the standard lightbulb process. 
Then for all n even, with a 1 as given in Q), 



sup 



X -n/2 
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H — for all n > 6, 
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where 



-r- 1 
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and for all n = 2m + 1 odd, 
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In the even case, as A n! 2. n decays exponentially last to zero, the variance a 2 is ol order n and the bound ([7]). 
therefore, ol order l/\/n; analogous remarks hold for the case where n is odd. 

We now more formally describe the lightbulb process on n bulbs with k stages. With n € N fixed and 
s = (si, . . . , Sk) with s r € {0, . . . , n} for r = 1, . . . , fc, we will let X s = {X^ : r = 0, 1, . . . , k, j = 1, . . . , n} 
denote a collection ol Bernoulli variables. The initial state olthe bulbs is given deterministically by {Xqj, j = 
1, . . . , n}, which will be taken to be state zero, that is, all bulbs off, unless specifically stated otherwise; non 
zero initial conditions are considered only in the Appendix. For r G {1, . . . , fc} the components ol the 'switch 
variables' X s have the interpretation that 



1 il the status ol bulb j is changed at stage r, 
otherwise. 



At stage r, s r ol the n bulbs are chosen uniformly to have their status changed, and the stages are 



independent ol each other. Hence, 
distribution ol X s is given by 



with e = {e r j}i< r <:k,i<j<n & n array ol {0,1} valued variables, the 



P(X S 



n;u(:)~ 





if yi 7 -— i &rj — s r, T — 1, . . . , k 



otherwise. 



(9) 



Clearly, the vectors of stage r switch variables (X r i, . . . , X rn ) are exchangeable and the marginal distribution 
ol the components X r j are Bernoulli with success probability s r /n. In general, for j = 1, . . . , n, the variables 



X; 



mod 2 and X= 



(10) 



are the indicator that bulb j is on at the terminal time fc, and the total number ol bulbs on at that time, 
respectively. For the standard lightbulb process we will write X and X for X n and X n respectively. 

The lightbulb process, where the individual states ol the n bulbs evolve according to the same marginal 
Markov chain, is a special case ol a class ol multivariate chains studied in QT] , known as composition Markov 
chains ol multinomial type. As shown in such chains admit explicit lull spectral decompositions, and in 
particular, the transition matrices for the stages ol the lightbulb process can be simultaneously diagonalized 
by a Hadamard matrix. These properties were put to use in [7] for the calculation of the moments needed to 
compute the mean and variance of X . Here we put these same properties to somewhat more arduous work, 
the calculation of moments of fourth order. 

That no higher order moments are required for the derivation ol a finite sample bound holding for all 
n is one distinct advantage ol the technique we apply here, Stein's method for the normal distribution, 
brought to life in the seminal monograph [TU]. By contrast, the method ol moments requires the calculation 
and appropriate convergence ol moments ol all orders, and obtains only convergence in distribution. Stein's 
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method for the normal is based on the characterization of the normal distribution in [S] , which states that 
Z is a standard normal variable if and only if 

E[Zg(Z)] = E[g'(Z)} (11) 

for all absolutely continuous functions g for which these expectations exist. The idea behind Stein's method 
is that if a mean zero, variance one random variable W is close in distribution to Z, then W will satisfy (|11[) 
approximately. Hence, to gauge the proximity of W to Z for a given test function h, one can evaluate the 
difference Eh(W) — Nh, where Nh = Eh(Z), by solving the Stein equation 

f(w) - wf{w) = h(w) - Nh 

for / and evaluating E[f'(W) — Wf(W)]. A priori it may appear that an evaluation of E[f'(W) — Wf(W)] 
would be more difficult than that for Eh(W) — Nh. However, the former form may be handled through 
couplings. 

Here we consider size bias couplings to evaluate E[f'(W) — Wf(W)]. Given a nonnegative random 
variable Y with positive finite mean fi = EY, we say Y s has the F-size bias distribution if P[Y S G dy] = 
(y/fi)P[Y € dy], or more formally, if 

E[Yg(Y)) = /j,E[g(Y s )] for all bounded continuous functions g. (12) 

The use of size bias couplings in Stein's method was introduced in [TJ, where it was applied to derive 
bounds of order cr -1 / 2 for the normal approximation to the number of local maxima Y of a random function 
on a graph, where a 2 = Var(Y). In [5] the method was extended to multivariate normal approximations, 
and the rate was improved to c -1 , for the expectation of smooth functions of a vector Y recording the 
number of edges with certain fixed degrees in a random graph. In [4] the method was used to give bounds in 
the Kolmogorov distance of order er -1 for various functions on graphs and permutations, and in [B] for two 
problems in the theory of coverage processes, with bounds of this same order. A more complete treatment 
of Stein's method and its applications can be found in |3]. 

Here we prove and apply Theorem l2.1l for bounded, monotone size bias couplings, which requires that the 
random variable Y of interest, and a random variable Y s , having the F-size bias distribution, be constructed 
on a common space such that for some nonnegative constant B, 

Y <Y S <Y + B 

with probability one. Loosely speaking, Theorem 12.11 says that given any such coupling of Y and Y s on a 
common space, an upper bound on the Kolmogorov distance between the distribution of Y and the normal 
can computed in terms of EY, Var(Y), B and the quantity 

A = yJVeLr(E(Y s - Y\Y)). (13) 

Theorem 12.11 is based on a concentration type inequality provided in Lemma 12.11 

For the standard lightbulb process, a size bias coupling of X to X s is achieved in the even case by the 
construction, for each i = 1, . . . , n, of a collection X 1 from the given X as follows. Recalling (fTDj) . where for 
s = n we have k = n, if Xi = 1, that is, if bulb i is on at the terminal time, we set X* = X. Otherwise, let 
J be uniformly chosen from all j for which X n / 2 j = 1 — A„/ 2 ,i an d let X* be the same as X but that the 
values of X n / 2>i and X n / 2 .j are interchanged. Let X % be the number of bulbs on at the terminal time when 
applying the switch variables X*. Then, with / uniformly chosen from 1, . . . ,n, the variable X s — X 1 has 
the X size bias distribution, essentially due to the fact, shown in Lemma |3.2[ that 

£(X l ) = C(X.\Xi = 1). 

Due to the parity issue, to handle the odd case when n = 2m + 1, we first construct a coupling of X 
to a more symmetric variable V. The variable V is constructed by randomizing stages m and m + 1 in the 
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switch variables that yield X. In particular at stage m we add an additional switch with probability 1/2 
and, independently, at stage m + 1 we remove an existing switch with probability 1/2. A size bias coupling 
of V to V s can be achieved as in the even case, thus yielding a bound to the normal for X. We remark that 
the size biased couplings developed here are used in [5] to show the distribution of X, in both the even and 
odd cases, obey concentration of measure inequalities. 

In Section [5] we present Theorem l2.ll which gives a bound to the normal when a bounded, monotone size 
biased coupling can be constructed for a given X. Our coupling construction and the proof of the bound for 
the even case of the lightbulb process is given in Section I57T1 Symmetrization, that is, the construction of V 
from X, coupling constructions for V, and the proof of the bound for the odd case are given in Section [XU 
Calculations of the bounds on the variance A in (|13l) require estimates on A n ,b,s in ©• These estimates, 
based on the work of |11| , yield the spectral decomposition of the underlying transition matrices of the chain 
and are given in Section |4j Some of the more detailed calculations have been relegated to the Appendix. 



2 Bounded Monotone Couplings 

Theorem 12.11 for bounded monotone size bias couplings depends on the following lemma, which is in some 
sense the size bias version of Lemma 2.1 of [§]. With Y having mean fj, and variance a 2 , both finite and 
positive, with some slight abuse of notation in the definition of W , we set 

W=?—^ and W s = ^—^. (14) 
a a 

Lemma 2.1 Let Y be a nonnegative random variable with mean fi and variance a 2 , both finite and positive, 
and let Y s be given on the same space as Y , having the Y -size bias distribution, and satisfying Y s > Y with 
probability one. Then with W and W s given in \1J$ , for any z £ R and a > 0, 



Proof. For fixed z £ R let 



-E{W S - W)ls W ,- W<a \li z<w<z+a } < a. 



—a w < z 
f(w) — ( w — z — a z < w < z + 2a 
a w > z + 2a. 



Then, using < a for all w £ R, Var(W / ) = 1 and the Cauchy Schwarz inequality to obtain the first 

inequality, followed by definition ([M)) and the size bias relation (fT2|) . we have 



a > E(Wf(W)) 



l -E { Y-»)f( Y -» 
a 

^E{f{W s )- f(W)) 
a 

„w a ~w 



rW —W 

= ^E / f'(W + t)dt 

a Jo 

rW'-W 

> ^E ±{0<t<a}l{z<W<z+a}f'{W + t)dt, 

G JO 



where in the final inequality we have used that W s > W and f'(w) > for all w £ R. Noting that 
f'iyV + t) = l{ z <w+t<z+2a}, an( l that < t < a and z < W < z + a imply z<W + t<z + 2a, we have 

l{0<t<a}l{z<W<z+a}f (W + *) = l{0<t<a} l{z<W<z+a} , 
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and therefore obtain 

,j f ws - w 

°~ Jo 

= (min(a, W s - W)l {z<w<z+a} ) 

> — E(W S - W)lf W s_ w<a x lf z <w<z+a} ! 
(7 

as claimed. □ 
Theorem 2.1 Lei y be a nonnegative random variable with mean fi and variance a 2 , both finite and positive, 
and let Y s be given on the same space as Y , with the Y -size bias distribution, satisfying Y < Y s < Y + B 
with probability one, for some positive constant B. Then with W and W s given in {1$ , we have 

sup \P(W <z)- P(Z < z)\ < 4A + 0.82—^ + 5, 

where 



A = ^Var{E{Y s -Y\Y)) and 5 = B/a. (15) 

Proof. For z G K arbitrary, let h(w) — 1{ W < Z } and let f(w) be the unique bounded solution to the Stein 
equation 

f'(w)-wf(w)=h(w)-Nh, (16) 

where Nh = Eh(Z). Substituting W into (|T6"|) and using definition (fT^) and the size bias relation (fT2"|) yields 

E(h(W)-Nh) 

= E(f'(W)-Wf(W)) 

= E(f'(W)-^(f(W s )-f(W))) 

I j-W s -W \ 

= E (f'(W)(l - ^(W s - W)) - £ / (f(W + t)-f(W))dt). (17) 



a a 



o 



As compiled in [3] , we have the following bounds on the solution / from Lemma 2 in Chapter II of [TU] , 

< f(w) < V27T/4 and \f'(w)\<l, (18) 
and also, as previously noted in [5], as a consequence of (|18|) and the mean value theorem, we obtain 

\(w + t)f(w +t)- wf{w)\ < (\w\ + V^7r/4)|t|. (19) 
Noting that EY S = EY 2 / n by (HD, we find 

^-^) = 4f^-^ = i. 

cr a 2 \ [i J 

Therefore, taking expectation by conditioning and then applying (|18|) and the Cauchy Schwarz inequality, 
we bound the first term in (JTTJ) as 



E {f'(W)E (l - ^(W s - W)\W^ < ^y/Va,r(E(W s - W\W)) = -^A. 



6 



To bound the remaining term of (fTTj) . using ([T5]). we have 
„w s -w 



(f(W + t)-f(W))dt 

r W s -W 



pVV — vv />W S —W 

- [(W + t)f(W + 1) - Wf(W)]dt + H / (l {w +t< z} ~ l{w<z})dt 

a Jo a Jo 



(20) 



Applying to the first term in (j2"0")l . and using < W s — W < 8 and EW 2 — 1, shows that the absolute 
value of the expectation of this term is bounded by 



E 



w s -w 



\W\ 



tdt 



2a 



E\{W S - Wf{\W\ + Z f L )\ < £^(1 + ^f 1 ) < 



S 2 fi 



Taking the expectation of the absolute value of the second term in (|2"0")) , we obtain 



»E 



W-W 



(l{W+t<z} — l{W<z})dt 



»E 



W'-W 



l{z-t<W<z}dt 



< ~E(W S - W)l {z _ s<w < z} 



again using that < W s — W < 5 with probability 1. Lemma I2TT1 with a = 5 and z replaced by z — 5 shows 
this term can be no more than 5. Since z G R was arbitrary the proof is complete. □ 



3 Normal approximation of X 

The next lemma shows that the size bias distribution of a sum may be achieved by taking certain mixtures. 
The result is a special case of Lemma 2.1 of [5], but we give a short direct proof to make the paper more 
self-contained. 

Lemma 3.1 Suppose X is a sum of nontrivial exchangeable Bernoulli variables X\, . . . ,X n , and that for 
i G {1, . . . , n} the variables X\ 1 . . . , X l n have joint distribution 

C(Xl...,X i n )=C(X 1 ,...,X n \X i = l). 

Then 

n 

has the X-size bias distribution X s , as does the mixture X 1 when I is a random index with values in 
{1, . . . , n}, independent of all other variables. 

Proof. For i G {1, . . . , n}, we first need to show that X i satisfies (H2J), that is, that E^E^X*)] = E[Xg(X)} 
holds for a given bounded continuous g. Now, for such g 

n n 

E[Xg(X)} = E \ X 39{X)] = £ P[X S = mgWlXj = 1]. 
i=i i=i 

As exchangeability implies that E[g(X)\Xj = 1] = E[g{X)\Xi = 1] for all j = 1, . . . , n, we have 
E[Xg(X)] = [J2 P [Xj = 1] I E\g{X)\Xi = 1] = E[X]E[g{X% 
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proving the first claim. The second claim now follows from 

n n n 

E 9 (x I ) = J2E[g(x I ),i = i) = Y, E ^( xI )\ I ^W = ^)^T, E 9( xt W = ^ 

i—l i—1 2—1 

n n 

= ]T Eg{X s )P{I = i) = Eg{X s ) £ P(I = i) = Eg{X s ). 
t=l »=1 

□ 

3.1 Even case 

In this subsection we provide the proof of Theorem II .11 for even n. We begin by describing a coupling of X, 
the total number of bulbs on at the terminal time n in the standard lightbulb process, to a variable X s with 
the X-size bias distribution. Throughout, we let U(S) denote the uniform distribution over a finite set S. 
Theorem 3.1 With n £ N even, let the collection of switch variables X = {X r j : r,j — 1, . . . ,n} and X 
satisfy (0) and U0\) , respectively, with s = n. For every i = 1, . . . ,n let X 1 be given from X as follows. If 
Xi = l thenX.* = X. Otherwise, with £(J*'|X) = U{j : X n / 2 ,j = l-X n / 2 ,i}, letX i = {X l rj : r, j = 1, . . . , n} 



Y" 1 



id to x i = YJLi X) wht 



X r j r 7^ n/2 
X n/2 ,j r = n/2, j £ {£, J 4 } 

Xn/2^ r = n/2,j = J i , 



Then, with I uniformly chosen from {1, . . . , n} and independent of all other variables, the mixture X 1 = X s 
has the X-size bias distribution and satisfies 

X s -X = 21 {Xl=QtXjl=0} and X<X s <X + 2. (21) 

We prove Theorem 13.11 making use of a preliminary lemma, and also of the fact that 

P{X 3 = 0) = P{X 3 = 1) = 1 forallj = l,...,n. (22) 

The equalities in (|22|) follow from EX — n/2, itself implied by ([2} and that Ai.„ jn = 0, as noted earlier. 

Lemma 3.2 For all i = 1, . . . ,n, the collection of random variables X 1 constructed from X as specified in 
Theorem lff.il satisfies 

£(X J ) = C(X\Xi = 1). 

Proof. For a given i £ {1, . . . , n} let e = {e r j : r, j = 1, . . . , n} with e r j € {0, 1} for r, j — 1, . . . , n. First note 
that since X 1 = X when Xi = 1 we have 

P(X l = e) = P(X t = l)P(X l = e\X t = 1) + P(X* = 0)P(X* = e\X, = 0) 
= P{Xi = 1)P(X = e\X t = 1) + P(X t = 0)P(X* = e\Xi = 0), 

so the desired conclusion is equivalent to 

P(X* = e\Xi = 0) = P(X = e|Xi = 1). (23) 
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As the construction of X 1 preserves the number of switches in each stage r we may assume e r j = r 
for all r, as otherwise both sides of (|23l) are zero. If 53 e ™ = 0mod2 then the left hand side of (|23|) is zero 
since X\ — 1 by construction; similarly, the right hand side is zero as X = e implies Xj = 0. Hence we need 
only verify (|23[) assuming 

n n 

e r j = r for all r = 1, . . . , n, and e r i = 1 mod 2. (24) 

3=1 r=l 

Writing J for J 1 for simplicity, and letting e 4J denote the array e with coordinates e n /2 t i and & n /2,j 
interchanged, by (|22|) we have 

P(X 4 = e\Xi = 0) = 2P(X 4 = e,Xi = 0) 

- 2^P(X i = ^ = 0, J = j) 

n 

= 2^P(X = e^,X i = 0,J = i). 
i=i 

Note that when e n /2,i = e n/2,ji or equivalently, e^ 2 i = e^y 2 ^ then 

P(X = e^.JQ = 0, J = j) = P(X = e« J = j) (25) 

as both sides are zero, since on J — j we have X n / 2: i 7^ X n / 2 .j- Otherwise, e„y 2 ,i 7^ e n/2,ji an d by the second 
equality in (|2"I|). 

n n 

e rf = e ™ e «/2 i = e r i + 1 — e n /2,i = e r i + 1 = 0, with equalities modulo 2, 

r=l r=l 



so (|25|) holds again. Hence 

P(X l = e\Xi = 0) 

n 

= 2^P(X = e^',J = j) 
3=1 

n 

= 2^P(J = j|X = e lJ ')P(X = e 4J ) 



3=1 



—In n / \ 1 



s=l v 7 3=1 

where in the second to last equality we have used that e'' J ', as e, satisfies the first equality of (PMl) . and 
the distribution of X given by ©, and in the last equality that the sum of probabilities of the conditional 
distribution of J given an X configuration that satisfies the first equality of (|24[) for r = n/2 must sum to 
one. Now, again using the second equality in (j24|) . 

2 J] y s j = 2p (* = e ) = 2p ( x = e > x > = 1) = p ( x = e l x > = i). 

proving (|2U)) . and the lemma. □ 
We now present the proof of Theorem 13.11 
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Proof. That X s has the X-size bias distribution follows from Lemmas 13. II and 1 3.21 To prove the first equality 
in (|2Tj) , note that if X] — 1 then X 7 = X, hence in this case X s = X. Otherwise Xj — and the collection 
X 7 is constructed from X by interchanging the stage n/2, unequal, switch variables X n / 2 j and X n i 2j ji. 
If Xji = 1 then after the interchange X\ = 1 and Xj t — 0, yielding X s — X. If Xji = then after 
the interchange X\ = 1 and X^ = 1, yielding X s = X + 2. The second claim in (|21l) is an immediate 
consequence of the first. □ 
The following lemma shows that for the case at hand the variance of the conditional expectation term 
(|15p in Theorem 12.11 may be expressed in terms of quantities of the form 

9a!p,B ~ = ■ ■ • = X a+ p = 0, Xn = ■ ■ ■ = Xi a = 0, Xi, a+ i = ■ ■ ■ = Xi >a+ = 1), (26) 

the probability that when applying the switch pattern s = (si,...,Sfc) on n bulbs over k stages, bulbs 
numbered 1 though a + /3 terminate in the off position, and in some stage I € {1, . . . , k}, bulbs 1 through a 
receive switch variable 0, and bulbs numbered a + 1 through a + j3 receive switch variable 1. We keep the 
number of bulbs n implicit in the notation, and also suppress s when s = n, writing more simply g^\- 

Using the spectral decomposition in Section 0] to handle the probabilities in (|26p , we now provide an 
upper bound to the term l|15p when applying Theorem 12.11 for n even. For a given s = (si, . . . , Sfe) and 
l€{l,...,k} let 

S; = (8i,..., 8l-i,Sl+i,..., 8k), (27) 

the vector s with its I th component deleted. For ( ^ j we similarly let sij denote s with its I th and j th 
components deleted. 

Lemma 3.3 Let n be even and X and X s given by Theorem \3.1\ Then for n > 6, 

-=-111 



A < A where A = J Var(E(X s - X\X)) and A = — = + — + -e~" /2 . 

2y/n 2n 3 

Proof. We apply the construction of X s , and the conclusions, of Theorem 13.11 For notational simplicity let 
J 1 = J, so in particular, from (|21[) we have X s — X = 21{Xi=o,Xj=Q}- Expanding the indicator over the 
possible values of / and J, and then over the values of the switch variables X n /2,i and X n / 2 j yields 

n 

1{Xj=0,Xj=0} = J! 1 {X l =0,X j =0} 1 {I=i,J=j} 

n n 
= ^2 1 {X l =0.X j =0,X n/2ti =0} 1 {I=i,J=j} + 1 {X l =0,X j =0.X n /2,i = l} 1 {I=i,J=j} 

= ^2 1 {X t =0,X j =0,X n/2it =0,X n/2J = l}'i-{I= l ,J=j} + ^2 1 {X i =0,X j =0,X n/2it = l,X n/2d =0} 1 {I= l ,J=j} 

= 22]l {x . = o,x j =o,x„ /J , j =o,x„ /2ij =i}l{j=i,j=j} 1 

where the second to last equality holds almost surely, as the probability of the event {/ = i, J = j} is zero 
whenever X„/ 2 ,i and X n /2j agree, and the last inequality holds as the final expression is the sum of two 
terms which can be seen to be equal by reversing the roles of i and j. 

To obtain a tractable bound on the required variance we apply the inequality 

Vax(E(X s - X\X)) < X'av(E(X s - X\T)), (28) 

which holds when T is any c-algebra with respect to which X is measurable (see [5], for example). Here 
we let T be the er-algebra generated by X, the collection of all switch variables. The first indicator in the 
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final sum above, l{Xi=o,Xj=o,Jf„ /2 i= o.x n/2 j= i}i is measurable with respect to J 7 . For the second indicator, 
conditioning on T yields 

e (i {I=ii j =j} \r) =p(i = i,j = = jLi {Xn/2 ^ Xn/2 . } , 

as for any i, chosen with probability 1/n, there are n/2 choices for j satisfying the condition in the indicator. 
Hence, recalling from (l2~Tj) that X s — X = 21{Xj=o,x 7 =o}; we have 



E X s -X 



J 7 ) = U n where U n = \ ^ l{x,=o,x :J =o,x„ /2il =o,x„ /2 , J =i}, (29) 

and < Var(C/„) by ([25]). 

Taking the expectation of ?7 n in (|29p . by the exchangeability of the (71)2 terms in the sum, applying 
Corollary 14. II and recalling the notation defined in (|27|) . we obtain 



EU n = ^(n)29[T } = \ (1 - An, 2 ,n n/2 ) • (30) 

Squaring (|29[) in order to obtain the second moment of U n , we obtain a sum over indices with 
{h,h} n {ji, j 2 } = 0, so |{ii,i 2 , ji, J2H € {2,3,4}, and we may write 

^ = ^,2 + ^,3 + ^,4 where (31) 

U n,p = ^4 l{X il =0,X 31 =0,X„ /2 , il =0,X n/2 , :fl =l} 1 {X <2 =0,X i2 =0,X„ /2 ,i 3 =0,X n/2 , 3 . 2 =l}- 

{ii,»2}n{ji,i2}=0 

Beginning the calculation with the main term U% 4 , where all four indices are distinct, taking expectation 
using exchangeability, Corollary 14.11 yields 

EU "a = ^W4<& /2) = (^) (! - 2A„, 2 ,„„ /2 + A„, 4 ,„„ /2 ) . (32) 

With the inequalities over the summation in (|3"Tj) in force, the event \{ix, = 3 can only occur 

when 

a) ii ^ 12, ji = 32 or b) iy = i 2 , ji ^ j%- 
Now continuing to use Corollary 14 . 1 1 without further notice, case a) leads to a contribution of 

16 f . ( n /2) _ n-2 , . 

( n )392,l — , 2 V 1 + A «,l,n„ /2 - A n,2,n„ /2 ~ ^n,3,n n/2 ) , 



while in the same manner the contribution from case b) is 



16, , (n/2) _ n- 2 , , 

— \ n )59l,2 — ^ 2 \ _ A n,l,n n/2 — A„,2,ii n/2 + A ni 3,n„ /2 J . 

Totaling we find 

^n,3 = ^(l-A„, 2 ,„ n/2 ). (33) 

With the inequalities over the summation in (|3Tj) in force, the event \{ii, i 2 , ji, J 2 }| = 2 can only occur 
when i\ = 12 and j\ = 32- Hence, 

EU* = ^(n) 2 ^/ 2) = i (1 - A„, 2 , n „ /2 ) . (34) 
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Summing ([52]), (f3"31 and (JM]) we obtain 

= (~4^~) (- 1 ~ 2Xn ' 2 > n »/2 + A «:4,n n/2 ) + — (l ~ A„, 2 ,n„ /2 ) ■ 

Subtracting the square of the first moment, given in (|30)) . yields 

1 ( 2 \ 1 1 

Var(C7 n ) = — f 1 - - J (l - 2A„, 2 ,n„ /2 + A„, 4 ,„„ /2 ) + — (l - A„ i2 ,„„ /a ) - — (l - A„, 2 , n „ /2 ) 2 

= Yg (A n ,4,n„ /2 - K,2,n n/2 ) + ~J^T i 1 ~ 2A n , 2 ,n„ /2 + A„,4,n n/2 ) + ^ i 1 ~ An,2,n n/3 ) 

1 1 1 / 2 \ 1 1-n 

~ An + 4n 2 + 16 V V4 ' n "/ 2 V2,n„ /2 J 2n2 A„, 2 ,„„ /2 + ^ A„, 4 ,„„ /2 . 

Now applying Lemma 14.51 from Section |4] for n > 6 we obtain 

1 1 „ / 1 1 1 

^ 1- + + e \T7 + ~2 + F~ 

4n 4n z \ lb n z 8n 

The inequality \/a + b + c < \fa + Vb + y^c, holding for all nonnegative a, b and c, now yields the claim of 
the lemma. □ 

With all ingredients at hand, we may now prove the bound for even n. 
Proof of Theorem ] 1. 11 even case. The size biased coupling given in Theorem 13.11 satisfies the hypotheses of 
Theorem 12.11 with B = 2, by the second inequality in (l2~Tj) . Hence the result for the even case follows by 
applying Theorem 12. II with /i = n/2,6 = 2/ a and the bound Ao on Ao given in Lemma T3. 31 □ 

3.2 Odd case 

Now we move to the case where n — 2m + 1 is odd. Instead of directly forming a size biased coupling to X, 
we first couple X closely to a more symmetrical random variable V for which a coupling like the one in the 
even case may be applied. The variable V is constructed by randomizing stages m and m + 1. In particular 
at stage m wc add an additional switch with probability 1/2 and, independently at stage m + 1 we remove 
an existing switch with probability 1/2. 

Formally, let X = {X r j : r,j = 1, . . . , n} be a collection of switch variables with distribution given by (j9]) 
with s = n, and let X — X n be given by (JTUJ) with k = n. Let 

C{B m \X) = U{j : X mj = 0} and £(S m+1 |X) = U{j : X rn+lj = 1}, 

and let C m and C m+ i be symmetric Bernoulli variables, independent of X and of B m and B m+ i. Now let a 
collection of switch variables V = {V r j, r, j = 1, . . . ,n} be defined by 

X r j r ^ {m, m + 1} 

X mj r = m,j^B m 

C m r = m,j = B m (35) 

X m+ ij r = m + l,j'/ B m+ i 

C m+ i r = m+ l,j = B m+ i, 



V rj 



and set 



V = ^ Vi where v i = \Y1 Vr i mod 2 - ( 36 ) 

3=1 Vr=l / 
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In other words, in all stages other than m and m + 1 the switch variables that produce V are the ones from 
the given collection X. In stage m, the switch variable for all bulbs but bulb B m , chosen uniformly over all 
bulbs in that stage that were not toggled, are the ones given by X. The switch variable for B m in stage 
m, however, is set to C m , which takes the values and 1 equally likely. Hence with probability 1/2, one 
additional bulb in stage m is toggled. Similarly, in stage m + 1, the switch variable of bulb B m+ \, uniformly 
selected from all the bulbs that were toggled in that stage, is no longer toggled with probability 1/2. Since 
X and V differ in at most two switches, we have 

\X-V\<2. (37) 

We now show some basic facts about the distribution of the switch variables V, that as in even case 
EV — n/2, and verify the variance formula ([6|). In the following, let 

n(a, b) — (1, . . . , m — 1, m + a, m + b, m + 2, . . . , n) for a, b G {0, 1}. (38) 

Lemma 3.4 The collections of variables {V r j,j = 1, . . . , n} are mutually independent for r = 1, . . . ,n. For 
r 6 {m, m + 1}, 

C(V r i, . . . , V rn ) = — £(X m i, . . . , X mn ) + — £(X m +l,l, • ■ • , -Xm+l,n)j (39) 

and furthermore, with Vj given by A36}) , we have P(Vj = 1) = 1/2 for all j = l,...,n, EV — n/2 and 
Var(V) = cry as given in (0|). 

Proof. The first claim follows by the independence of {X r j,j = 1, ... ,71} for r = 1, . . . , n, and that {V r j,j = 
1, . . . , n} = {X r j, j — 1, . . . , n} for r ^ {m, m + 1}, and otherwise {V r j,j = 1, . . . , n} is given by randomizing 
{X r j, j = 1, . . . , n} independently of the remaining stages, and of stage 2m + 1 — r. 

To prove (|39|) . first consider the case r = m. Let ei, . . . , e„ € {0, 1}. Since (V m i, ■ • ■ , V mn ) = (X m i, . . . , X mn ) 
when C m — 0, and since C m is independent of {X m \, . . . , X mn ), we have 

P(V m i = ei, . . . , V mn = e„) 

= -^P(V m i = ei, . . . , V mn = e n \C m = 0) + -P(V ml = e 1: . . . , V mn = e„|C m = 1) 

= -^P{X m i = ei, . . . ,X mn = e n \C m = 0) + --P(Kni = ei ; ■ • ■ , Knn = e n \C m = 1) 

= -P(X, n i — ei, . . . ,I m „ = e„) + -P{V m i = ei, . . . , V^„ = e n \C m = 1). 

Hence (|3"9"|) is equivalent to 

-P(Kji = ei, . . . , F„ m = e n |C m = 1) = P(X m+ i,i = ei, . . . ,X m+1) „ = e n ). (40) 

Now assuming . ej = m + 1, as both sides of (|4TH) are zero otherwise, using the independence of C m from 
.B m and X we have 

P(V m i = ei, . . . , V^ m = e n \C m = 1) 

n 

= P{ y^ = ei, • • • , V mn = e„|B m = j, C m = l)P(S ro = j|C m = 1) 
3=1 

= Yl P ( X ^ =^ l +3, X mj = 0\B m = j, C m = l)P(B m = j) 

= P ( X ^ = e l' l ^j' = 0|# m = j)P(B m = 3) 

= P ( B " = 3\ X ™i = e h l ^3\ X m3 = 0)P{X ml =e h l^j, X mj = 0). 
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Now, since B m is chosen uniformly from the n — m variables taking the value in stage m, recalling 
J2j e j = m + 1, we have 

P(V m l = ei, . . . , V mn = e n \C m = 1) = ^ P(X m l = C[ , I ^ j, X m j = 0) 



n — m 

j:e j= l 



m + 1 / n \ I n 



n — m\mj \m + 1 



P(Xm+l,l — e l; • ■ • ) Xm+l,n ~ e n), 



proving (|4U1) . The argument for r = rn + 1 is essentially the same. 

Recalling the notation in by the distribution of V is the equal mixture of £(X S ) over the four 
cases s = n(a, b) for {a, 6} € {0, 1}. Hence, from ([2]), 

EV ~~Q X! (! ~ A «4,n(a,6)" 



2' 

r,t£{m,m+l} 



as 



7 ^n,l,n(o,6) 7 ( H ^n,l,s ] (A^ )lim + 2A n ,i, m A„,l, m +l + A^j : „ l+ i) 



n,l,m H~ -^n.l.m+l )/2)' = o, 

s^{m,m+l} 



since 



A 2m \ A 2(m + l)\ 

An l m + A n ,l,m+1 — 1 + ( 1 — 0. 

V « / V n J 

Clearly EV = n/2 implies P(V} = 1) = 1/2 

Next, by the conditional variance formula and ([2]) we obtain 

^ / / TL \ ^\ I f Ti Tl^ \ 

Vtu(V) = - J2 Var(X n(Q , )) + i?(X n(a , b) --j =- ^ ( -(1-A n , 2 ,„(a,b)) + — A n ,2,n(a,6)) , 

a,6G{0,l} a,be{0,l} 

and now applying (1/4) Y, a ,be{o,i} ^n,i,n(o,6) = A„,i jn , and a similar identity for A„ )2 ,n( a ,fe), yields the vari- 
ance cry in ([5]). □ 

We now present a size bias coupling for V. As in the even case, the variable V s is obtained by first 
constructing, for each i = 1, . . . , n, switch variables V 4 that satisfy 

£(V l ) =C(V\Vi = 1). (41) 

For a given i = 1, . . . , n, to construct V J , one first determines if Vi = 1. If so, set V = V. Otherwise, let M 
be a variable that chooses from the stages to and m + 1 uniformly and independently of V. As in this case 
Vi — 0, one may achieve V? = 1 by changing the switch variable Vui to 1 — Vm%- The coupling accomplishes 
this change in one of two possible ways. 

To introduce the first way, called a flip, we say a configuration e of binary switch variable values is feasible 
if P(V = e) ^ 0, that is, when 

n n 

e r j = r for r ^ {to, to + 1} and e r j € {m, m + 1} for r G {m, m + 1}. 
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If flipping Vmj to 1 — Vmi results in feasible configuration, then the flip is made with probability l/(m + 1). 
In other words, given e, r £ {m, m + 1} and i = 1, . . . , n let e" be the configuration with entries 

e s ; s 7^ r or I ^ i 
1 — e r i s = r and / = i. 

Defining V '* in like manner, the distribution of F l , the indicator that Vm« is flipped, is given by 

P(F l = 11 V, M) = — — 1(V M < J is feasible). (42) 
m + 1 

In the flip is unsuccessful, that is, if F l = 0, we perform an 'interchange' in stage M, much like the 
coupling in the even case. For a configuration e, r £ {to, to + 1} and i, j £ {1, . . . , n}, let 

e s i s ^ r or I ^ i 
e r i s = r,l£'{i,j} 

&ri & — I j 
&<p j S — f ' j / — t ^ 

that is, e 1 "'*'- 7 is the configuration e with the variables in the r, i and r, j positions interchanged. Now let J l 
be a random index with distribution given by 

£( J l \V, M) = U{j : V Mj + V Mi }- (43) 

Defining V ' l ' 3 in like manner, when F" 1 — 0, we interchange Vm% with Vm j*- 
Hence, overall the configuration V z is specified by 

V t = l 

V t = 0, F* = 1 (44) 
M,i,r y. = Qj pi = q 

The following theorem shows that V 1 satisfies (|41l) . 

Before stating the theorem, we claim that if e is feasible, then for r £ {m, m+1}, the configuration e r ' 1 is 
feasible if and only if X^yj e rj = m - Both e and e r ' 1 are feasible if and only if e rj + e ri £ {m, m + 1} and 
Ylj^ti erj+1— G r i £ {to,to+1}. The claim now follows noting that {to— e r i, m+1— e r i}n{TO—l+e r i, m+e r i} = 
{to — 1, m} H {m, m + 1} = {to}. Summarizing, for a feasible configuration e, 

e r,J is feasible if and only if e r j — m that is, if and only if e r i = Y^j—x e rj — m. (45) 

In the following we denote F l and J 1 by F and J respectively for simplicity. 
Theorem 3.2 Let V be constructed from X as in \35)) , let M be a random variable uniformly distributed 
over {m, m + 1} and independent V, and for i £ {1, . . . , n}, given V and M let F and J have distributions 
as specified in J^l?| ) and respectively. Then V 1 given by J^[ ) satisfies J^i[ j. 

Further, letting V 1 = Ylj—x V* where 




\r=l 



w/ien / is uniformly chosen from {1, . . . , n} independent of all other variables, the mixture V 1 — V s has the 
V -size bias distribution and satisfies 

V s - V = l {Vr=0 ,F=i} + 21 m=0 ,x. r =o,F=o} and V < < V + 2. (46) 
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Proof. As in the even case, since V* = V when Vi = 1 and P(Vi = 1) = 1/2, the proof that V 1 has the 
V-size bias distribution follows by verifying, for all i = 1, . . . , n, that 

P(V l = e\Vi = 0) = P(V = e|y 4 = 1) or equivalent^ P(V 4 = e, V = 0) = P(V = e, V = 1) (47) 

holds for all feasible configurations e satisfying ^"=1 e ™ = 1 niod2. 

Recalling that M, F and J are the stage selected, the indicator that the bit in question in the V config- 
uration is flipped, and the index with which an interchange is to be performed, respectively, we find 

P(V* =e,V = 0) 

= P(V* = e, Vi = 0, F = 1) + P(V* = e, V- = 0, F = 0) 

■m+l / n 

= ^] P(V j = e, V 7 ; = 0, M = r, F = 1) + ^ P(V J = e, V = 0, J = j, M = r,F = 0) 



3=1 



m+l 



5^ P(V = e r '\ Vi = 0, M = r,P = 1) + ^ P(V = e r, ' J , V» = 0, J = j, M = r,F = 0) 



3=1 



m+l 



= ( P ( V = ^ M = r » F = !) + S = e " J ' J = J '' M = r ' F = °) I ( 48 ) 

3=1 

where in the final equality we have used that X)r=i e ™ = lm°d2. 

Let X)j=i e i-3 = m + b, where necessarily b £ {0, 1}. For the first term in the r th summand of (|48|) . by (|45[) . 
(jUJ), the independence of M from V and that P(V = e) has the same value for any feasible configuration, 
we have 

P(V = e r,i ,M = r,F = 1) 

= P(P = 1|V = e", M = r)P(V = e r <\ M = r) = -P(V = e) 1 ^" = b \ (49) 

2 m + l 



Next, for the second term in the r th summand of (|4"5)) . we have 

n 

P(V = e r ' id , J = j,M = r,F = 0) 

3 = 1 

n 

= ^P(J = i|P = 0, V = o' ' ■'. .\/ = r)P(F = 0|V = (•" •'. M = r)P(V = e rAi ,M = r) 

3 = 1 

1 ™ 

= -P(V = e)^P(J =i|P = 0,V = e r,i,J ',M = r)P(P = 0|V = e r ' lJ ,Af = r), (50) 

3=1 

the final equality holding by again using that M is uniformly chosen over {m,m + 1}, independent of all 
other variables. 

To handle the summand terms in (|50p . note that the index J chooses uniformly over the m + b(l — e^;/'"') + 
(1 — b)e^ indices whose stage r switch value is opposite to e r r ^' J , that is, 

P(J = j\F = 0, V = e^,M = r) = -l^^p - . . = - -. (51) 



+ 6(1 - elf J ) + (1 - ft)e™ m + be rl + (1 - &)(1 - e ri ) 



Again using e^ lJ = e rj -, by (gSJ and (|4"2"j). 



■ ■ ?77 

P(P = 01 V = e r ' hl ,M = r) = l(e H = 1 - 6) + l(e H = 6). (52) 

m+l 
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Taking the product of ([51"]) and (j52")l yields 



m + 6e„ + (1 - 6)(1 - e ri ) \^ = 1 " &) + ^TT 1( ^ = 6) J = i^+l 1 ^ * e ^ (53) 
Hence may be written as 

P(V = e) -A 1 / m 



•>, „, +17 51 x ( e « ^ = 2 P(V = e) ( ^TT 1(e " = 6) + 1(e " = 1 ~ b) 

j=l 



Adding this factor to (1491) we find that the r summand of (|48j) equals 

1 P(V = e) 1(e " = &) + JP(V = e) f-^rlCert = b) + l(e ri = 1 - b)) = Jp(V = e). 
2 ?7i + 12 \m+ 1 / 2 

Hence, adding up the two summands r — m and r = m + 1 in ()48j) yields (|47p . as 

P(V J = e, Vi = 0) = P(V = e) = P(V = e, Vi = 1), 

where we have applied X)"=i e " = 1 m °d 2 for the final equality. 

Lastly, to prove (|4l)|) . first note that by construction V s — V when Vi = 1. When Vi — and P = 
then y s equals either V or V + 2, depending on whether Aj takes the value 1 or 0, respectively, as in the 
even case. When Vj = and F = 1 the status of bulb J is changed from off to on and the status of all other 
bulbs remain unchanged, hence V s — V = 1. □ 

With the coupling of V to V s now in hand, we prove a bound to the normal for V. We recall from 
Lemma [3"H that EV = n/2 and Var(V) = cr^, given in (JBJ. 
Theorem 3.3 Ifn= 2m + 1, an odd number, and V is given by I136\) , then 



sup 



^2 77, 2 

< — -5-A1 + 1.64^- + — for all n>7, 
2cr^ a v 



where o~y and Ai are given in ^) and {3J) respectively. 

Proof. We apply Theorem 12.11 to the coupling construction given in Theorem 13.21 As S = 2 /ay from (1461) , 
and EV — n/2 from Lemma 13.41 it remains only to show that yVar (E (V s — V\V)) < Ai for all n > 7. 
As in the even case, we obtain an upper bound by conditioning on the cr-algebra T generated by V, with 
respect to which V is measurable. Again for notational simplicity we drop the superscripts on F and J 
unless clarity demands them. 

Taking conditional expectation in (l46|) yields 

E(V S - V\T) = Cn + 2Cn where ( n = E (l {Vl=Q , F=1} |P) and £ n = E (l {Vl=0>Vj = 0jF=0} \F) . (54) 

Letting 

A = Var(C„) and P = Var(£„), 
since |Cov(£„,£ n )| < V 1 AB by the Cauchy-Schwarz inequality yields, we have 

Var(P(F s - V\F)) <A + 4VAB + AB. (55) 



First computing a bound on A, expanding the indicator defining in (J54J) over the possible values of M 
and I yields 



m+l n 



l{Vr=0,P'=l} — /] l{M=r,f=i,Vj=0,F=l}- 
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Recall that conditional on V, configurations for which V" are feasible, that is, by PS")) , those for which 
Vrj — m i ar e flipped with probability l/(m + 1). Now, since / and M are chosen uniformly and 
independently over indices {1, . . . , n} and {m, m+1} respectively, taking conditional expectation with respect 
to T we obtain 



-y m+l n ( j 

— V V 1 { Vi = 0, V V rj = m > 

2ra(m + l) ^ ^ I ^ 3 | 

r—m i—1 I j^i I 



(56) 



To determine E( n we now demonstrate that the events {Vi = 0} and {Ej^; ^rj = ™} are independent 
for r £ {m, m + 1}. Note that for each fixed i, 

WXXi = TOj = PlVi = o,J2v rj = m ] +p v; = i,^v>,= 

so it suffices to prove that 

P[Vi = Q,J2 V rj= m ) =P \ Vi = l,J2 V rj= m ) ( 57 ) 
V jfr J \ j#« / 

since P(Vi = 0) = 1/2, by Lemma 41 However, as the map that sends e to e r ' 1 is a bijection between the 
two events in (|57|) , since all feasible configurations have equal probability the equality in (IBTl) holds, and 
{Vi = 0} and {Ej^i Vrj — m} are therefore independent. 
Now, using that 

F E^ =m = P [^Vrj = m,C r = \ + P l^2v rj = m,C r = 1 



= lp(v ri = o\c r = o) + ±p(y ri = i\c r = i) 



c r = 1 



l/m + l m + 1 \ m+l 



(58) 



2 \ n n J n 

and that P(Vi — 0) — 1/2, taking expectation and using independence in (|56p we have 

1 m+l n , 

U 2n(m + l)^^ 2n 2n K ' 

r—m i—1 

Turning to the second moment of £ n , squaring the right hand side of (1561) we obtain 



I n n f 



4n 2 (m 

We decompose An 2 (m + l) 2 Cn m t° the following four terms: 

1. For r = t, i = I we obtain CL = E"l™ EIU = 0, E;^ ^ = m} 

2. For r = t, * ^ Z we obtain Cl 2 = E"l + m E*i = 0, V, = 0, £ j¥i = m, £\y V>, = m} 
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3. For r ^ t, i = I we obtain <^, 3 = Er/t E?=i 1 W = °> Ej# = TO > Ej#j = m > 

4. For r ^ t, t + I we obtain £ 4 = E r#i E^z 1{^< = 0, V, = 0, ^ V rj = m, £ J?4J - ™} 
For the expectation of Cn.i: using independence as in the calculation of E( n and ([55)1 . we obtain 

771+1 71 

^Cn,i =EE P W = °> E ^ = m) = m + 1. (60) 

r—m i—1 j^i 

Turning to Q 2 we h& v e 

{V< = 0, F = 0, V rj = m, ^ V rj = m} 

= (J {^ = 0,^=0, ^ = m-b, V ri =b,V rl = &}, 
&e{o,i} jg{i,t} 

and hence may write 

_ m+1 f ] 

C",2= E Cn,2, 6 where C',2,6 = E E 1 i** = °' V * = °> E V rj =m-b,V ri = b,V r i = b\. (61) 

6e{0,l} r=m i^i [ jg{i,2} J 

Hence, to compute 2 b we sum 

P I Vi = 0, Vi = 0, 53 V rj = m - b, V ri = b, V H = b 

\ mi,t) 

= P Ivi = 0, Vi = 0, ^ ^ = m - 6, Kni - 6, y mi = 6 

= P ^ = 0,^ = 0,5^^ =m + b, V mi = b,V m i = b 

= P(Vi = 0, Vi = 0, V mi = b, V m i = b,C m = b) 

= 51 P(Vi = 0,V l = 0,V mi = b,V m i = b,C m = b,C m+1 = a) 

ae{0,l} 

= \ E P(^i = 0,F 2 = 0,y ml = 6 ! F m2 = 6|(7 ro = 6 ! (7 ro+1 = a) ! 

a£{0,l} 

where in the second line we used the exchangeability of stages m and m + 1 in V, and for the final equality 
that i I, and exchangeability again. 

Recalling definition (13"5|) . for a, 6, G {0, 1} easily we have that 

£(V|C m = b, C m+ i = a) = £(X n ( fei0 )), 

and so, applying definition (I26[) . we obtain 

P (14 = 0, V 2 = 0, T/ ml = b, V m2 = 6|C m = 6, C m+1 = a) - 4-2&,26,n(&,a)- 

Hence, recalling (|6TJ) , 

p>2 _ Wa M too'* 

■^71,2 ~ 2^ S , 2_26,26 1 n(6,o)- 

a,6e{0,l} 
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Now turning to E( 2 3 , by again considering the bijection that maps e to e" we may conclude that 

P = 0, v n = "*> ^2V tj =m\ = P F; = 1, ^ V rj = m, ^ V tj = 

and hence, by the independence of the switch variables in stages r 7^ t, as provided by Lemma l3.4[ we obtain 
(m + l) 2 



and therefore 



P ( ^ = = m = 2P Vi = 0, Vrj =m,^2 V tj = m 



<3 = ^^- (63) 

2 



Arguing similarly to compute £ 2 4 , the identity 
(to+ l) 2 



= p \Y. Vr 3 = m >J2 ^ = m = 4P K = 0, = 0, ^ V rj = to, J] 7y = 



yields 



2 _ (n) 2 (m + l) . , 

^C„, 4 ^ • ( 64 ) 



Combining (JBDJ, (JB2JI, (JM]), ||BD and O yields 

a = s(c 2 )-£;(c„) 2 

= 1 ffm I 1) I (n)2 V o (m) 1 ( m + 1 ) 2 ■ (nh(m + lA 1 

4n 2 (m+l) 2 I 2 ^ ^-ibflb^h,^ n 2n 2 \ 4n 2 ' 

To obtain a bound on A first note that for any 6 we have n(6, 0) m = n m+1 and n(6, l) m = n m , and that 
the m th component of n(6, a) is to + 6. Hence when 6 = Corollary 14. II yields 

(m) __ ( |(1 + 2A n ,i, nm+1 + A„,2,n m+ i) ^rl)? 2 a = 

#2,0,n(0,a) ~~ ] 1/1,9), 4- X ) ("+1)2 ' -1 



and when 6=1 that 



f(l — 2A ni i )Ilro+1 + A n , 2 ,n m+ i) < '7nV>'' 2 a — 



(m) jA n,l,n m +i ~r ^n,2,n m + i/ (71)2 

»0,2,n(l,o) _ 1 1/1 9X , a \ (m+l) 3 " „_i 

Summing these four terms, and using that (A n . 2 ,n m + A„,2,n m+ i)/2 = A, lj2 ,n,„, as defined in ([S]), yields 
,M _ 1 (a , o\ , 01 \( m + 1 ) 2 _/i,T ,(m + l) 2 



2^ 9 2 -2b,2b,n(b.a) = 7 ( 4 + ^ A ™,2,n m + 2A„, 2 ,n m+1 J — 7"7 - (1 + A„, 2 ,„„J 7-— 

o,6S{0,l} l j2 W2 



Hence, 
.4 



f 1 -1 \ j (")a A w + 1 )2 , t (m + l)2 \ (to + 1) 2 , (n) 2 (TO + l) 2 \ 1 
(m + 1) + -7— — --r h A„,2,n m — p— v + 



4n 2 (ro + l) 2 \ 2 \ (n) 2 ' ' m (71)2 / Ji 2 n 2 J An 2 

1 / 1 to I n— 1 \ to — 



4n 2 \to + 1 2(m + l) n 2n J 8n 2 (m+l)' 

3m + 2 to — 

8n 3 (TO+l) + 8n 2 (TO + l) ™' 2,nm ' 



20 



Applying Lemma Ff~51 we obtain 

3m + 2 _3_ 

~ 8n 3 (m + 1) + 8n 2 " 8n 3 + 8n 2 ' ( j 

Next, expanding the indicator in (|54|) to obtain B = Var(£ n ), and recalling that an interchange is 
performed when F = 0, we have 

1{Vj=0,Vj=0,F=0} 
m+l n 

= ^ ^ l{V 4 =0,V i =0} l{I=i,J=j,M=r,F=0} 
r=m i,j=l 
m+l n 

= ^ ^ (l{V i =0,V J =0,V ri =0} + l{V i =0,V r 3 =0,Vri=l}) ^-{I=i,J=j,M=r,F=0} 
r=m i,j=l 

m+l 

= (•'■{Vi=0 ) V i =0,V r i=0,Vri=l} + l{V5=0,V5=0,V ri =l,Vri=0}) l{/=i,J=j,AF=r,F=0} 

r— m z^j 

m+l 

— ^ ^ ^ l{v i =0,V r j=0,V r4 =0,V r 3=l}l{7=i,J=j,M=r,F=0}! 
r—m i^j 

where the last inequality holds due to the symmetry between i and j in the penultimate sum, as in the even 
case. 

As V is ^-measurable, recalling the definition £ n = J3(1{Vj=o,Vj=o,f=o} m (|54[) . taking conditional 
expectation yields 

m+l 

£« = 2 ^ ^ 1(^=0,^=0^^=0, K J = l}£'(l{/=i,J=j,M=r,F=0}l^ r )- (66) 

r—m i^j 

To compute the conditional expectation in the sum, using the independence of / from M and V, and M 
from V, we obtain 

£(1{ J=i, J=j,M=r,F=0} K ) 

= P( J = j|M = r, I = i, F = 0, V)P(F = 0\I =i,M = r, V)P(J = i\M = r, V)P(M = r|V) 

= — P(J = j|M =r,I = i,F = 0,V)P(F = 0\I = i,M = r,V) 
2n 

= ■z-PiJ* = i|Af = r,^ = Q,V)P(F i = 0\M = r.V) 
2n 

1 ljVri ± Vrj} 
In m + l 

the final equality since the product of ([52"]) and (|FT|) is given by (|53p . Hence, as the first indicator in the sum 
in (po) specihes that V>j = 0, V r j = 1, we obtain 

m+l 

61 = —7 — TT / , / , l{Vi=0, V5=0,Vr4=0, Vrf=l}- (67) 

n(m + l) ^— ' ' 1 

v 7 r=m 

Taking the expectation of £„, by the exchangeability in V of stages m and m + l, we have 

Efa) = , \ n E p{ y- = °> f j = °> y ™ = °> ^ = 

n(m + 1) ' 
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Exhausting over the possible values of C m = b and C m +i = a for a,b £ {0, 1}, we obtain 

p(Yi = o, Vj = o, v mi = o, v mj = i) 

- Y = °> *J = °> ^ = °> ^ = 1 |C m = 6, C ro+ i = a) 



4 

a,6G{0,l} 



1 \- (m) 

I. fl,l,n(6,o) 



4 

a,bG{0,l} 

Now, by using Corollary 14. 11 

P > n-1 M _ n — 1 , -r- v m(m + l) _ m , T x 

^" - 2fa+T) a J- i} 5 M,n( 6 ,a) - ^7" ^T) (1 - V2,nJ ^ " ^ (1 " A„,2,n m ) • (69) 

To calculate squaring in (|67|) we obtain a sum over indices r, t € {m, m + 1} and i%, i%^2\iH with 
H 7^ ill *2 7^ J2i and therefore i 2 , ji, = P fo r f> S {2, 3, 4}. Hence we may write 

£ = Y &,P Where &,p = Y €n,p,q and 

pe{i,2,3,4} <?e{i,2} 
&,p,q = n 2( m + i)2 Y Y Mvi 1 =oy Jl =oy r , il =o,v rdl =i}^v i2 =oy j2 =oy t , i2 =o,v t , :i2 =i}' 

\{r,t}\=q Ki».-S2,n,i2>!=P 

First consider the case p = 4, that is, {?i , «2 , ji , J2} | =4. We begin with ^ 41 , that is, when r — t. Using 
again that {C r — b, C n — r — a} partitions the space for a,b € {0, 1}, we have 

P{V h = 0, V l2 = 0, V h = 0, V j2 = 0, V rM = 0, V r>i2 = 0, V r , h = 1, V rd2 = 1, a = 6, C„_ r = a) 

= -P^ = 0, V l2 = 0, V h = 0, K A = 0, VrM = 0, K, l2 = 0, V r . n = 1, V r , j2 = l|C r = 6, C„-r = a) 

p(Vi = 0, v 2 = 0, u 3 = o, v 4 = 0, v rA = 0, = 0, k, 3 = 1, v;,4 = i|a = 6, C n _ r = a) 



4 
L 

^2,2,11(6, a) 



1 » 



Summing over the four distinct indices ii,«2,.7i,j2 and applying Corollary 14 . 1 1 we obtain 

(n) m+l 
E (C,4,,l) = 4 n 2( m+ l)2 Y Y #2,2,n(&,a) 
V ' a,6e{0,l} r=m 

2(^)4 /-. cyr , T >, (m) 2 (m + 1)2 

V ZA n,2,n,„ + -"Vi,4,n m J 



16n 2 (m + l) 2 ' ' m ' ' (71)4 

In order to proceed further we need to consider functions such as those in (|21)1) . but which specify 
switch variables in two stages, rather than in only one. In particular, for s = (si, . . . , Sfe), distinct stages 
r, t 6 {1, . . . , k} and bulbs i\ < ■ ■ ■ < i a and ji < ■ ■ ■ < jp in {1, . . . , n}, let 

\k,a 1 ),...,{i a ,a a );(hM),-,(Jp,bii)],s ( 71 ) 
= P(Xi 1 = 0, . . . , Xi Q = 0, Xj 1 — 0, . . . , = 0, I r ,ti = ai, . . . , X ri j Q = a a , X t j 1 = 61, . . . , X J/3 = &s) 

when X has switch pattern s. We note that it is not required that . . . ,i a } PI {ji, . . . ^ 0- The 
functions in (fTTj) give the probability that bulbs {ii, . . . , i a } U{j'i, . . . , jp} are off, and that the switch values 
for bulbs ii, . . . ,i a and ji, ■ ■ ■ ,jp assume the prescribed values in stages r and t, respectively. 
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Now moving to 4 2 , considering first r = m, i = m+l and arguing similarly as for the £ n ,4,i, 
P(V h = 0, y 42 = 0, V n = 0, V j2 = 0, V^, 4l = 0, V m+lti2 = 0, F mjl = f , V m+1 , ]2 = 1, C ro = 6, C m+1 = a) 

= = 0, ^ 2 = 0, ^ = 0, = 0, Kn,l = 0, Kn,2 = 1, W,3 = 0, V m+1A = l\C m = b, C m + 1 = a) 
1 (m,m+l) 

- 4%l,0),(2,l)i(3,0),(4,l)],n(6,a)' 

applying the definition (I7f I) . Summing over the four distinct indices, noting that the case r = m + l,t = m 
contributes equally, using (| 103[) from the Appendix we obtain 

fiYf 2 1 - ( n )4 \ " (m,m+l) 

^Un,4,2j - 2n 2 fm + f 1 ) 2 ^ 9 [(l>0),(2,l);(3,0),(4,l)], n (M) 



(n) 4 



a : be{0,l} 

2 



2n 2 (m + 1) : 



( (m (n)! )3 ) (( 2m - 1 ) 2 + 2 ( 2w - 1 )V2,n m , m+1 + A„,4, Ilm , m+1 ). (72) 



Summing ^(^4,1) and S(^ 4;2 ) from ([70)1 and (|72l respectively, and simplifying the latter using n = 2m + l 
yields 

jYi ( ffl — ■ 1 ^ 

E (&^ = Sn 2 (m + 1) ^ ~ 2A "- 2 - n '" + A "- 4 -"-) ( 73 ) 
+ 2n 2 (r+l) 2 (n) 4 ( (2W - 1)2 + 2 (2 ™ - 1} A - 2 — + 

m 2 (m- 1) i - - > (m) 2 (2m - 1) (m) 2 . (m)j 

8n 2 (m + l) v ' ' ' 8n 3 4n 3 ' "*' m+1 2n 2 (n) 4 ' ' + 

Moving to p — 3 and the evaluation of 3 (J ), under the constraints i\ ^ j\, i 2 7^ j 2 , sums over indices 
that satisfy 12, Ji , } | = 3 can be partitioned into the following four cases: 

= h,3i h, 2. i x ^ i 2 , ji = 32, 3. i x = j 2 , «2 ^ ji, 4. i 2 = ji, i\ 7^ j 2 , 

Subscripting by these subcases, here and similarly in what follows, we write 

£n,3,g = £ra,3,g,l + £n,3,g,2 + £n,3,g,3 + £n,3,g,4' 

Starting with q = 1, for i?£ 2 311) we have 

P(VS = 0, ^ = 0, V^ 2 = 0, V r ,t = 0, V rtjl = 1, V;, J2 = 1, C r = b, C n - r = a) 

= ~P(Vi = 0, V 2 = 0, V3 - 0, V r ,l = 0, V r ,2 = 1, K,3 = l|C r = 6, C n _ r = a) 

1 (r) 
~~ 4"l,2,n(6,a)' 



while for E£ 2 



n,3,l,2i 

PfVi! = 0, V i2 = 0, Vj = 0, K,m - 0, V tM = 0, = 1, C r = b, C n _ r = a) 
= ip(F a = 0, V 2 = 0, V 3 = 0, K,i = 0, K,2 = 0, K r , 3 = l\C r = b, C n _ r = a) 

1 (r) 

4"2,l,n(6,a)" 
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Under q = 1, cases 3 and 4 have probability zero, so summing and applying Corollary 14.11 yields 

(n) m+l 

£(£n,3,l) = E (Cn.3,1,1 + ^1,3,1,2) = ^ n 2( m T f\2 (^lj.nCM) + 32,l,n(b,a; 

V ' r=TOo,6e{0,l} 

_ ( n )3 ( ( m ) 1 (m) 

_ 2n 2 (m + l) 2 ^ V ffl .2,n(6,a) 92,l,n(b,a] 
y ' a,be{0,l} 



m(2m — 1) , — 
An z (m + 1) y 



(74) 



Moving to q = 2, for ££^321 we compute 



P(Fi = 0, V h = 0, = 0, V r ,i = 0, V t ,i = 0, 7 rJl = 1, Vt, h = 1, <? r = 6, C t = a) 

1 = 0, U 2 = 0, U 3 = 0, V r ,i = 0, V r , 2 = 1, H,i = 0, VJ.s = l|C r =b,C t = a) 



1 >,*) 



4 
1 

- 4-9[(l,6),(2,l);(l,0),(3,l)],n(6,a) 

while similarly, with factors of 1 /4 for 3 2 2 we have 



P(Vi = 0, V 2 = 0, K 3 = 0, V r ,\ = 1, K r , 2 = 0, Vt,i = 1, Vt, 3 = 0|a = 6, C t = a) = gf^^ 



(2,0);(l,l),(3,0)],n(6,a)' 



for ^n,3,2,3 we have 



(r,t) 

[(l,0),(2,l);(l,l),(3,0)],n(6,a)' 



P{Vi = 0,V 2 = 0, U 3 = 0, V r ,i = 0, U r , 2 = 1, Vi,i = 1, VJ, 3 = 0|C r = b, C t = a) = g 
and for E£ 2 3 2 4 we consider 

P(U X = 0, U 2 = 0, V 3 = 0, K,i = 1, V r ,2 = 0, Vt,i = 0, Vt, 3 = l|a = 6, C t = a) = J4m),(2,o);(i,o),(3,i)],„(M)- 
Summing up these q = 2 case terms and applying (|108[) in the Appendix, we obtain 

^(£71,3,2) — £(£n,3,2,l + £n,3,2,2 + £n,3,2,3 + £n,3,2,4) 

(n)a 



4n 2 (m + l) 2 

v 7 a,be{0,l} r^i 



EV- / (r,t) (r,t) 
^[(l,0),(2,l);(l,0),(3,l)],n(6,a) + 9'[(l,l),(2,0) ; (l,l),(3,0)],n(6,o) 



1 (>",*) 1 ( r >*) ^ 

+ %L,0),(2,l);(l,l),(3,0)],n(6,a) + %l,l),(2,0) ; (l,0),(3,l)],n(6,a)J 

(n) 3 /4(m + 1)| (2m - 1) 2 \ m(2m - 1) 



4n 2 (m + l) 2 V {n)\ ) 2n z 



(75) 



Now summing (|75l) and (1741 we obtain 

„ i0 m(2m — 1) / ^- s m(2m — I) 



It remains to consider E£~ l2 - Under the constraints ii 7^ ji,ia ^ j'2, sums over indices that satisfy 
*2, Ji, J2H = 2 can be partitioned into the following two cases: 

1- h =ii,h =h and 2. i x = j 2 , ji = i 2 . 
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Consider first q=l. Under Case 1, by consideration of the term P(Vi = 0, V 3 ■ = 0, V r j — 0, V r ,j = l,C r = 
b, Ct — a) and the application of Corollary 14. 11 we obtain 

in) m+l ^ 

^n,2,l,l = 4n 2( m + 1 )2 J2 <i,n(6,a) = 2 n 2 (m + 1) ^ ~ ' (7?) 

Case 2 has zero probability. 

Now let q — 2. For stages r t, under Case 1 we consider 

P(Vi = 0, Vj = 0, = 0, V r .,j = 1, Vf^ = 0, V t} j = l,C r = b, C t = a) = ^3[( r i*o) : (2,i);(i:0),(2,i)]:n(6,a)- 

Now, summing over the two ways one can achieve r ^ t, which yield identical terms, and using (|113|) from 
the Appendix, we obtain 

Fp i _ (nh y- (r,t) 



2n 2 (m + l) 2 ^ y [(l,0),(2,l);(l,0),(2,l)],n(f ) ,a) 
V ; a,f,G{0,l} 



1 + 2A ni i jnm m + 1 + A rl! 2,n mj 



2n 2 (m + l) 2 v ' "'■"-' n —+i 1 '"M.-».™+w (n) 2 

^— g (1 + 2A„ i i ! n m _ m+1 + A n! 2,n m , m+ i)- (78) 



Under Case 2 we consider 



P( Vl = 0, ^ = 0, V r ,i = 0, 7 nj = 1, U M = 1, V t j =0,C r = b, C t = a) = ^ [ ( 53 ) , (a , 1) . (1 , 1 ), (2)0)]>n(6ia) 
and similarly arrive at 

^£n,2,2,2 = 4^3 ~ 2 ^n,l,n m , m+ i + K,2,n m , m+1 )' (79) 

Summing up the terms ([77]). (1751 and (JTSJ) yields 
Finally, recalling 

B = Var(^) - £ (£ 4 + £ 3 + ^,2) - (^n) 2 , 
collecting terms (|75|) . (|75)l . and (j8"0)l and simplifying we obtain 

m(8m + 3) m(m 2 + m — 1)— (771 + 1)2. 

m 2 (m — 1) — ("1)2 \ m 2 — 2 

+ 8n 2 (m + l) A " Anm + 2r^j; ^ 4 ' nm ' m+1 ~ Ir^"' 2 '"™' 

Applying Lemma T4. 51 we find for n>7 

„ m(8m + 3) / 1 1 1 1 \ „ m(8m + 3) e -n 
B < - 5 + — + — + — + tt~ ~ e~ n < 



8n 3 V 16 16n 64 64n 2 J 8n 3 11 
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Now, using (frj5"j). 



< 



< 



3 m(8m + 3) /m(8m + 3)e-" / 3 e~ n e~ n e 



n 3 8n 3 V Sn 3 8n 2 V 8n 3 11 V 8n 2 11 



5m + 1 / 1 e~™ / 3 e 



8n 3 V An 8n 2 V 8n 3 11 V 88n 



3 — n 



< 



5m + 1 + ./ J_ e -»/a + J— e""/ 2 + — 
8n 3 V 32n 3 V 88?i 3 8n 

5m + 1 e~" /2 e"" 



2n 3 / 2 8n ' 
Hence, from (|55[) . again using n > 7, 

Var(E(V s - V|.F)) < A + 45 + 4a/AB 

, 3 e"™\ /m(8m + 3) e" n \ /5m + 1 e~"/ 2 



e 



8n 3 8n 2 J \ 8n 3 11 J \ 8n 3 2n 3 / 2 8n 

1 e""/ 2 



< 

~ n 8 

Taking square root and using y/a + b < ^fa + \fb now yields the upper bound © , thus completing the proof. 
□ 

We now provide a bound for the normal approximation of X in the odd case. We remark that using V, 
fewer error terms, and therefore a smaller bound, results when standardizing X as in Theorem ll.il that is, 
not by its own mean and variance but by the (exponentially close) mean and variance of the closely coupled 
V. 

Proof of Theorem\TJ$ odd n. Letting W = (X - n/2)/a v and W v = (V - n/2)/a v , recalling \X - V\ < 2 
from (f3"T|) . we have 

\W-W V \ = \X-V\/a v <2/a v . (81) 

With = P(Z < z) and 



^ n -r , „ . n 2 
C v = — g-Ai + 1.64— + — 

by (|8ip and Theorem 13.31 we obtain 



P(VF < Z) - $0) < P(V7y - 2/o~y <Z)- $(z) 

= P(W V <z + 2/a v ) - $(z + 2/cry) + $(z + 2/cry) - $(z) 
< Cy + 2/(cryV / 27T). 

As a corresponding lower bound can be similarly demonstrated, the claim is shown. □ 



4 Spectral Decomposition 

In |11| the lightbulb chain was analyzed as a composition chain of multinomial type. Such chains in general 
are based on a d x d Markov transition matrix P that describes the transition of a single particle in a 
system of n identical particles, a subset of which is selected uniformly to undergo transition at each time 
step according to P. 
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In the case of the lightbulb chain there are d = 2 states and the transition matrix P of a single bulb is 
given by 



P = 



1 

1 



where we let eo = (1,0)' and e\ = (0, 1)' denote the and 1 states of the bulb, for off and on, respectively. 
With b 6 {0, 1, . . . , n} let P n ,b,s be the 2 b x 2 b transition matrix of a subset of size b of the n total lightbulbs 
when s of the n bulbs are selected uniformly to be switched. Letting P n o a — 1 f° r ah n and s, and J 2 the 
2x2 identity matrix, for n > 1 the matrix P n fi, s is given recursively by 

P n .b,s = - (P ® P„_i,6_i, s -i) + (1 - -) (J 2 ® Pn-i,6-l,») for 6 e {1, . . . ,n}, 

as any particular bulb among the b in the subset considered is selected with probability s/n to undergo 
transition according to P, leaving the s — 1 remaining switches to be distributed over the remaining 6—1 
of n — 1 bulbs, and with probability 1 — s/n the bulb is left unchanged, leaving all the s switches to be 
distributed. 

The transition matrix P is easily diagonalizable by the orthogonal matrix T as 



P = T'TT where T = 



1 

V2 



1 1 

-1 1 



hence P n ,b,s is diagonalized by 



where r nj o,s = 1, and is given the recursion 



i h T'T r . 



and T = 



1 T 



1 

-1 



r n , 6>s = - (r <g> r„_ M _ M _i) + (1 - -) (i 2 <g> r„_i >6 _i jS ) for b e {1, . . . , n}. 



(82) 



(83) 



(84) 



The next result describes the diagonal matrices T ni b, s more explicitly in terms of a sequence of vectors of 
length 2 b for all 6 > 1 defined through the recursion 



a b = (a b _i,a b _i + l b _i) for b > 2, with ai = (0, 1), 
where 1& = (1, . . . , 1) is of size 2 b . For example, 

ai = (0,l), a 2 = (0,1,1,2) and as = (0,1,1,2, 1,2,2,3). 



(85) 



(86) 



Letting a n be the n term of the vector a b for any b satisfying 2 > n results in a well defined sequence 
ai,a 2 , ■ ■ .. 

Lemma 4.1 For n € {0, 1, ...,}, b, s G {0, . . . , n} 7 and X n ,b,s given by (Qp, the matrix T n ^, s in WW satisfies 

^n.b,s — diag( K X n ^ ai )S , . . . , \n,a 2 b ,s)- 

In particular, with 2 b-i the vector of zeros of length 2 b ~ 1 , for b > 1 



r n ,6,s — diag(X n7aiiS , . . . , A n! a 2() _ 1 lS , 2 b-i ) + diag(0 2 b-i , A„ !ai +i, 
For instance, from (f86| , for b — 2 we have 

r n ,2,s = diag(A„^o,s, A n! i, s , A ni i iS , A n>2iS ), 

and for 6 = 3, 



i A n)0 6-1 +!,«)• 



fn,3,s — diag(A n! o,s, A raj i iS , A ni i !S , A n! 2, s , A n! i^ s , A n! 2,s, A ni2)S , A U) 3 ]S ). 
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Proof. As a\ = we have r„.o, s = 1 = A„.o, s , so the lemma is true for 6 = 0. For the inductive step, assuming 
the lemma is true for b — 1, by (|84|) and the definition of T from (|82j). it suffices to verify 



— A n _i ,a,s— 1 + (1 )^n— l.a.s = A„, Q , S and 

n n 

s s 

' — A n — 1 ,o s— 1 + (1 )A„-1 a s — A n a+l,S) 



for all a = 0, 1, . . .. To prove the first equality, by ([T]) we have 

s \ , ft s n _ s ( a \t ^t( s ~ l )t , /-, S ^/«V Nt (s)t 



-An-!..-! + (1 - -)A„_ 1>0>S = -E U ("^HV + (!--) E U (- 2 ) r , 
n n n T^,\ t ) l n ~ l )t n fr£\tj {n — l)t 



t=o 



£/ \n(n-l)t V n/(n-l) t 



t=o 



*/ V (») 



t+1 



E ( a ) (-2)' ^ ^ s ~ 1 + 11 ~ ^ 



t=o 

' a 



t=o 



t f (*)*(»-*) 



t) V (») 



t=0 

A n _ a _ s . 

The second equality can be shown in similar, though slightly more involved, fashion. □ 
We note that expresses these eigenvalues in terms of the hypergeometric function. 
If the k stages of the process 1, . . . , k use switches s = (si, . . . , Sfc), then since the matrices P n ,b,s, s £ 

{0, 1, . . . , n} are simultaneously diagonalizable by (|83[) . the transition matrix P n ,b,s for any subset of b bulbs 

can be diagonalized as 

k 

P n ,b,s = II Pn *^ = ® h T'T nAs ® b T = ® b T'diag(A„, Ql , s , . . . , A„ >Q26 , S ) <g> 6 T, (87) 
i=i 

where A n! a !S is given in ([T]) and 

If 7T is a permutation of {1, . . . , k} let 7r(s) = (s^m, • • • > s n(k))- As all matrices r nj b,s are diagonal, from 
we have r nj & )S = r„ & 7r ( s ), and now from (|87|l we have the following result. 

Lemma 4.2 TTie distribution of the lightbulb chain is independent of the order in which the switch variables 
s are applied, that is, for all permutations ir, 

Pn,b,s — Pn,b,-rv(s)- 

The following lemma helps us compute probabilities such as g^ B s in (f26]). For j £ {0, 1, . . .} let fibj be 
the 2 b x 2 b diagonal matrix in the variables Xk, k £ {0, 1, . . .} given by 

Q b ,j = diag(x ai +j, • ■ • , x a +j ), and set u b = ® b Tef b and w b = ® b Tef b , (89) 
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where we recall eo = (1, 0)' and e\ = (0, 1)'. Note that for 6 = 1 we have 

til = Teo = — /=(1, —1)' and u>i = Tei = — 1)'. 
V 2 v 2 



Lemma 4.3 Lei i € {0, 1, . . .} and Q t — &t,o> an d suppose that for some vector !J t £ I 2 , that 



3=0 



(90) 



(91) 



holds for t = b — 1 with some sequence a(j), j = 0, . . . , b — 1. TTien ZioZrfs /or £ = 6 w/ien replacing v t by 
Vb = u\®Vt and a(j) by 



b-J 



*(?")+ (f)otf-l). 



and /or i = 6 lu/ien replacing vt by Vb — ui\ ® v t and a(j) by 



a w (j) 



b-J 



a(j)-(J)a(j-l). 



Proof. By (185|) we may write 



and by (|§0|) we have 



^6-1,0 
Qf,-!,! 



Hence, when Ufc = <g> «fc_i = (u^-i! — 1, {,-i)'/v / 2j we obtain 



v' b flbUb 



2 (Vb-l^b-l.O^fc-l + Wfc-l^h-l,l w b-l) 



2 <> 



1 

1 6-1 
3=0 
1 6 

3=0 
1 6 

3=0 

b 



6-1 
J 

6-1 
3 

6-1 



1 6-1 

3=0 
b 



6-1 

3 



a(j>i+i 



3 = 1 

6-1 

i-i 



a(j - 1) hz. 



a(j) + ( | Ki - 1) ) *j 



( J°«C?>j. 



(92) 



(93) 



as claimed. The proof is essentially the same, using (|90]) . when Uf, = iui ® Ub_i = (^_ 1 ,«[_ 1 )'/v2. □ 
In Corollary 14 . 1 1 we express the functions g^\ g in (|26[) using the corresponding conditional probabilities 

/a l/J)B = P(X t = 0, t = 1, . . . , a + p\X ,i = 0, i = 1, . . . , a, X 0>i = 1, * = a + 1, . . . , a + (3) (94) 



that are the subject of the next lemma. As with g^a s , we drop the dependence on s when s — n 
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Lemma 4.4 For given a, ft > 0, setting b — a + ft the probability f a ,8,s i n is given by 

1 b fb\ 

fa,0,s = ^ ( ; ) a <x,p(j)*n,3,B, (95) 



2 b 



where a a ,o(J) — 1 /or aZZ a > and 



Oo^O') = ( —1T~) a a,p-i(j) - ( 7 ) a a,/3-i(j - 1) for all a > 0, ft > 1. (96) 



6 J-""-^' \b 

Proof. Using exchangeability for the first equality, extracting the relevant component of the fc-step transition 
matrix and applying (|87[) we obtain 



fa,?,. = ((e' 1 ) 8/3 ® (e' )® a ) P„, fc ,.ef = ((e[)^ ® (e' )® a ) ® fc T'T n>b , s ® b Tef b = ^r„ )M u 6 , 

where u a and u>^ are given as in (1891) and Vb — wp ® u a . Hence, with fib — flb t o as in (|89p. the result follows 
from 

1 6 /o\ 

v'b^bu b = 7^^2 I . W/jGO^i- (97) 



2 b 
j'=o 



We first prove the case /3 = by induction in a; note in this case Vb = Ub- Equality (1971) holds with 
a a,o(j) = 1 for a = 0, as both sides equal x in this case. Assuming that (|97|) holds for some a > with 
a a,o(j) = 1) then (l9"2l of Lemma T4.3I implies that (|9T|) holds for a + 1 and ft — with 

Oa+i,o(j) = aQ >°G0 + act >°(j ~ X ) = L 

Hence (|9"?| holds for all a > and /3 = with a at o(j) = 1. Similarly, assuming now that (|9"T)) holds for 
«a,/3-i(i) with nonnegative a,ft—l, we have that (|97|) holds with a at p(j) given by (|96|) by (|93|) of Lemma 
14.31 thus completing the induction. □ 
As our computations involve only moments up to fourth order, we highlight these particular special cases 
of Lemma 14.41 in the following Corollary. 

Corollary 4.1 For a,ft>0, the probability g^K s in A26\) is given by 



1a B a = fa,P,a,Pa,p[B,l) where p at g{S,l) = — . 

(n)a+9 



(98) 



and f a ,/3, s is given by \95)) . For < a + ft < 4, the sequences in \95\) specialize to 

oo,o(j) = 1 

oo,i GO = (-1) J a i,o(i) = 1 

Oo.aO') = (^ 1 ) 3 > «i,i(j) = 1 - J o^d a 2,o(j) = 1 

ao,aC?) - ai, 2 0') = (-1) 3 "(1- 2j/3), a a ,i = 1 - 2j/3 and a 3>0 (j) = 1 

and 

o , 4 (i) = 01,3(7) = (-l) 3 '(l-j/2), a 2>2 = l-i(4-j)/3, a 3 ,i = l-j/2 and a 4 , (i) = 1. 
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Proof. By Lemma EL~2l that is, the fact that the switch variables can be applied in any order, conditioning 
on the values of the switch variables in stage I yields the same result as assuming these values as initial 
conditions in stage 0, and applying the switch pattern s;, that is, s skipping stage /. Hence the first claim 
in follows, as the first factor is the probability of the given event conditioned on the values in stage I, 
while the second factor is the probability of the conditioning event, as 

p(x ll = --- = x la = o,x l , a+1 = --- = x ltEt+ p = i) = [J ( - - j - - ) Yl ( — j— — ) =p aif3 (s,l). 



i=0 v 7 i=0 



The specific forms of the sequences a a ,p{j) for < a + f3 < 4 follow directly from the initial condition 
and recursion in Lemma 14.41 □ 
Applying Corollary 14. 11 we obtain, for example, the formulas 

(0 1 m , \ \ \ , ( n - sihsi 

52,1,8 — q\ L + A 'i,l,si — A n ,2,Sj — -V3,sJ , 

8 (nj 3 

and 

(0 1 n o\ j_ \ ^ (s;) 2 (n-s;) 2 

52,2,s — 16 ZA ».2, Si + ^n,4,sj ■ 

Lastly we present the bounds on products of eigenvalues of the chain used to handle the variance term 
(fT5)l when applying Theorem 1 1.1 1 to the lightbulb chain. 

Lemma 4.5 For all even n > 6, and s = n n «, or for all odd n = 2m + 1 > 7 and s £ {n m , n m _|_i, n miTn ,-|_i} ; 

|A„ )2 ,s| < e~™ and |A„, 4 , S | < ^e~™. 
Moreover, with X n ,b,s as given in f^), /or aZZ odd n = 2m + 1 > 7 we have 

|A„,2,n,J < e~ n and |A„,4,nJ < \^T r ' \ 

Proof. The following calculations slightly generalize the arguments of [7]. Let m be given by n = 2m, and 
n = 2m + 1, in the even and odd cases, respectively. For n > 2 consider the second degree polynomial (cf. 
©) 

/ 2 (x) = l-- + i% 0<x<n. 
n (n)2 

It is simple to verify that f2{x) achieves its global minimum value of — l/(n— 1) at n/2, and that f?,(x) has 
exactly two roots, at (n + \/n)/2 and (n — y/n)/2. Hence, as f$(x) < for all x between these roots, and 
additionally, as (x — l)/(n — 1) < x/n for all x € [0, n], we obtain the bound 

j _J_ for x E [ ""^ "+v^ ] 
1/2^)1 < | (1 ^ )2 forx£[ ^^ jen[0)4 (99) 

For £ 6 K let |_^J an d \x] denote the greatest integer less than or equal to x, and the smallest integer 
greater than or equal to x. In both the even and odd cases let 

t = i \ g 1 ' * " ' L 2 * J ^ ' ' m ' 

so that 

iA„, 2 ,n m , m+1 i = [ n 1/2(^)1 ] ( n i/ 2 ( S )i ) n i/^i 
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If either of the roots (n — y/n)/2 or (n + \fn)/2 is an integer then equality holds as both expressions above 
are zero. Now assuming neither value is an integer, the product is over disjoint indices. 

Applying the bound [n/2 — x\ + \n/2 + x\ = n and 1 — 2(n — s)/n = —(1 — 2s/n) yields 



' * n (i-!) 2 (n^ n 




1 



x 2 \ 

2s x - 



where |t| is the cardinality oft. 

Using 1 — x < e~ x and that [x\ > x — 1 on the first product, we obtain the bound 

\\ n2n + J < ^ e -f(L IL i^J)(L IL ^J+l)/2'j 4 e -|t|log(n-l) < e -(n-2Vn-l+2/Vn+|t|log(n-l))_ 

To control |t|, note that as fx] < x + 1, we have 

iti = L^j-r^i-i>^-3. 

As log 63 > 4, the result follows for all n > 64 from 

-2\^-l + 2/Vn+|t|log(n-l) > -2-y/n - 1 + 4(V« - 3) = 2^ - 13 > 0, 

and by direct verification for all odd 7 < n < 63, thus proving the claimed inequality for X n ,2.n m m+1 for all 
odd n > 7. 

By (|99]l we have |/2(x)| < 1 for x £ [0,n], and hence 

\K,2,n m I = \h(m + l)A n>2 ,„ mim+1 | < | A„, 2 ,n m , m+ i | , (100) 

so the claimed bound holds for A n ,2,n m > an d likewise for \ n ,2,n m+1 , for all odd n > 7. The claim for 
A n .2,n m = (A n! 2.n m + Ki,2.n m+1 ) /2 now follows. Direct verification of the claim for A ni 2,n m for all even 
6 < n < 62 now also completes the proof for all cases involving A„.2, s - 

Now turning to A n ,4 jS , for all n > 6 consider the fourth degree polynomial 

/4( ,) = l-^ + ^-«i + ^)i, 0<x<n. 
n (n) 2 (n) 3 (n) 4 

It can be checked that the four roots of f±{x) are given by 

n ± VV2\/3n 2 - 9n + 8 + 3n - 4 , n ± V-^VSn 2 - 9n + 8 + 3n - 4 
.ti± = and x 2 ± = , 

and that additionally the three roots to the cubic equation f[(x) = occur at 



yi = n/2 and y 2 ± — (n ± \fin — 4)/2. 

These roots satisfy 

< x X - < V2- < x 2 - < vi < x 2 + < J/2+ < xi+ < n. 
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To obtain a bound over the interval [xi_,xi+], evaluating fi{x) at its critical values we obtain 

3 3 
(n - l)(n-3) ~ (n- 3) 2 



/4(yi) = — tyt^ ^ < — ^ and 



2(3n 2 - 9n + 8) 6 
Ml/2±J - ~ n ( n _ i)( n _ 2 )(n-3) " ~(n-3) 2 ' 

To bound /4(x) by / 2 (x) in the remaining part of [0,n], write 

2 / \ / 16x(n — x)p(x) 



/ 2 (x) - / 4 (x) 



(n — l) 2 n 2 (n 2 — 5ra + 6) 
where 

p(x) = (An - 6)x 2 + (6n - 4n 2 )x + n 3 - 2n 2 + n. 
The roots of the quadratic p(x) are given by 



n 1 n(n — 2) 
Z±= 2 ± 2V^3- 
As 5n 2 — 15n + 12 > for all n, we have (2n — 3)(3n — 4) > n(n — 2), and therefore 

(y2\f?,n 2 - 9n + 8 + 3n - 4~) (2n - 3) > n(n - 2). 

Dividing by 2n ~ 3 and taking square roots demonstrates that 

X]_ < z_ < z + < xi + . 
Hence p(x) is nonnegative on the complement of z + ], and we obtain 

lf 4 (x)|<l torxe ^' Xl +] (101) 

1/41 jl " I forx^x^xx+LxG^n]. [WL) 

Now write for short 



C(n) = \J V^^Sn 2 - 9n + 8 + 3n - 4 so that Xi_ = - — 

Using ([TUT]) . (IM1) and that |s - n/2\ > y/n/2 whenever s < (n - C(n))/2, 1 - x < e~ x for x > and finally 
[xj > x — 1, we obtain 



1^1 L^^J /L^^j ' ' 



n iAn,4, s i< n / 2 2 (^)< n f 1 --) < e -<^<tfc»>+i>+c<»>v-H*c(»v»). (102) 

Now let 

f r n — C(n)-, . n + C(n) . } , , 
u= H ••• »L ^\>\{^,m + l}. 

Using A n! 4^ s = \ n ,4,n-s, (| 1 1 [) . |u| > C(n) — 3 and (|102[) . we have that 

|A„,4,n mi ,„ + i| 

L 2 J n / L 2 J \ / g \ C(n) — 3 



\ 2 

A n.4,s 



n iAn,4 )S i n i a ».4, s i n n 

< e -2(ri-2(C(n) + l)+C(n) 2 /n+2C(n)/n) + (C(n)-3)(log(6)-21og(n-3)) 

_ g-n . e -(n-4(C(n) + l)+2C(n) 2 /™+4C(n)/n+(C(n)-3)(21og(n-3)-log(6))) 



(n - 3) 2 
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Since for n > 96 we have 

n - 4(C(n) + 1) > log(2) and (C(n) - 3) (2 log(n - 3) - log(6)) > 0, 

we conclude that |A„.4 !nm m+1 | < e~"/2 for all such n > 96, and direct verification shows this same bound 
holds for all 7 < n < 95. Arguing as in (|100j) . the claimed inequality therefore holds for \ n .4,n m and A„.4 !nm+1 , 
and hence also for A n! 4. nm = (A n! 4. nm + \ n .4,n m+1 )/2 for all odd n > 7, as well as for X n ,4.n n/2 for all even 
n > 8. The proof is completed by directly verifying the bound for \ n .4,-n. n/2 f° r 1 = 6. D 

Appendix 

We work out some of the detailed calculations in the proof of Theorem 13.31 First, we show 

E (m,m+1) 
ff [(l.< 



,0),(2,l);(3,0),(4,l)] 1 n(6,o) 
a,b£{0,l} 

2 



( (m (n)4 )3 ) " + " 1 ) A "' 2 ' n —+^ + A ».4,n m , m+1 ) 



(103) 



To handle this sum, first write 

E(m,m+1) 
%,0),(2,l);(3,0),(4,l)],n(6,o) 

a,6e{0,l} 

= J]] J]] -P (^i =0,i= 1, . . . , 4, X m i = 0, X m 2 = 1, = ei, X m 4 = e2 

a,bG{0,l} ei,e 2 ,e 3 ,e 4 e{0,l} 

= e 37^m+l,2 = e A,^m+l,3 = 0, X m+1]4 = 1) . (104) 

We may write the probability in the sum as the conditional probability of {Xi = 0,i = 1,...,4} given 
{(X mll X m2 ,X m3 ,X m4 ) = (0, l,ei,e 2 ), (X m+1<1 , X m+1<2 , X m+1>3 , X m+1A ) = (e 3 , e 4 , 0, 1)}, multiplied by the 
unconditional probability of this last event. By Lemma 14^1 that is, that the switch variables may be applied 
in any order, conditioning the first event on the values of the switch variables for these four bulbs in stages m 
and m+ 1 is equivalent to conditioning on the combined event {(Xoi, X02, X03, Xq^) — (e^, €4 + 1, e±, ei + 1)} 
in, say, an initial stage 0, and running the lightbulb process in the remaining n — 2 stages, with, recalling 
our extension of the notation defined in (|27| . pattern n m , m+ i — (1, . . . , m — 1, m + 2, . . . , n). 
For the conditional probability, by (|94|) we may write 



P {Xi = 0, i = 1, . . . , A\(X 01 ,X 02 ,X 03 ,X 04 ) = (e 3 , e 4 + 1, e x , e 2 + 1)) = / Q ^,n m , m+1 (105) 
where a + j3 — 4 and 

a = l(e 3 = 0) + l(e 4 + 1 = 0) + l(ei = 0) + l(e 2 + 1 = 0). 

Having conditioned on the switch variables in stages m and m+ 1, we note that this conditional probability 
does not depend on a or b. 

For the unconditional probability, we have 

P {{X m i, X m 2, X m3 , X m 4) — (0, 1, ei, e 2 ), (X m+ i^i, X m+ i : 2, ^m+i,3) -^m+1,4) = ( e 3,e4,0,l)) 

\m+ 1 - 6) 3 _ ei _ e2 (m + fe)i +ei+C2 \ / (m + 1 - a) 3 _ e3 _ e4 (m + a)i +e3+e4 



(n) 4 / \ (n)4 

as the stages are independent, and the event requires 1 + e\ + e 2 switches to be drawn from m + & in stage 
m, and 1 + e 3 + e4 switches to be drawn from m + a in stage m + 1. Hence, changing the order of summation 
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in (|104j) . and recalling that the conditional probability f a ,p,n m m+1 does not depend on a, b, we have 



(m,m+l) X n 

o) — / Ja 



,0),(2,l);(3,0),(4,l)],n(6,a) ~ J a,fl,n mim + 1 W a , 
a,fee{0,l} a=0 



where 

= J2 J2 

E a ,f> {a,b}G{0,l} 



(m + 1 - 6)3_ ei -e 2 ( TO + &)i+ei+e 2 \ / (m + 1 - a) 3 _e 3 -e 4 (to + a) 



(n) 4 



Il+e3+ei 



(n) 4 



with 



Letting 



E a ,p = {(e 3 , e 4 , e x , e 2 ) : l(e 3 = 0) + l(e 4 + 1 = 0) + l(e x = 0) + l(e 2 + 1 = 0) = a}. 



P*(d)= E 



6e{o,i} 



(to + 1 - b) k -d(m + b)d 
{n)k 



we may write 



«J a ,i3 = E P 4 ^ 1 + ei + e 2)^4(l + e 3 + e 4 ). 



(106) 



(107) 



Writing out the sets E a ,[} required, we have that 

£ 4 ,o = {(0,1,0,1)}, 

E 3 ,x = {(0, 1, 0, 0), (0, 1, 1, 1), (0, 0, 0, 1), (1, 1, 0, 1)}, 

E 2 ,2 = {(0,0, 0,0), (0,0, 1,1), (0,1, 1,0), (1,0, 0,1), (1,1, 0,0), (1,1, 1,1)}, 

£1,3 = {(0,0,1,0), (1,0,0,0), (1,0,1,1), (1,1,1,0)} and E 0A = {(1, 0, 1, 0)}. 

As Pk(d) — Pk(k — d), and since Ep^ a is obtained from E a ^ by replacing each entry of each vector by its 
binary complement, 

w a ,/3 = 22 P4(3 -e\ - e 2 )p 4 (3 - e 3 - e 4 ) = 2J M 1 + ei + e 2 )p 4 (l + e 3 + e 4 ) = ius.q.. 



Now we calculate the required weights. As the one element of E^o satisfies (ei + e 2 , e 3 + e 4 ) = (1, 1), we 
have 



Wifi = p 4 (2) 2 



2(to) 2 (to + 1)2 



(n) 4 



Next, as the four elements in E% t x yield the vectors (1, 0), (1, 2), (0, 1), (2, 1) for (ei + e 2 , e 3 + e 4 ), we have 



tfi 8l i = 2p 4 (2)p 4 (l) + 2p 4 (2)p 4 (3) = 4p 4 (2)p 4 (l) = 8 



(to) 2 (to + 1)2 
(n) 4 



(to + l) 3 m + (to) 3 (to + 1) 
(n) 4 
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Lastly, as the six vectors in £2,2 yield (0,0), (0,2), (1, 1), (1, 1), (2,0), (0,0) for (ei + e2,&3 + e 4 ), we have 
w 2 ,2 = 2p 4 (l) 2 + 2p 4 (l)p 4 (3) + 2 J34 (2) 2 = 4p 4 (l) 2 + 2p 4 (2) 2 



(to + 1)3771 + (771)3(771 + 1) 
(n) 4 

Now using w a ^p — wp >a , applying Corollary 14. II in (|106[) we obtain 

fc,m+l) 

,0),(2,l);(3,0),(4,l)],n(6,o) 



(to) 2 (to + 1)2 



(n) 4 



E(m,m+1) 
%,C 



a,6e{0,l} 

— (1 + 6A„,2,n mim+ i + A ni 4 inrn>m+ i) 



(to) 2 (to + 1)2 

(77)4 



+ (1 — A n) 4 inm m+1 ) 



(777)2(771 + 1)2 \ / (777 + 1)3777 + (777)3(777 + 1) 



(n) 4 



+ Yg(l ~~ 2A ni 2,n m , m+ i + A n) 4 1 n TOjm+1 ) I 4 



Simplification yields equation (11031) . 
Next we demonstrate that 



(n) 4 

(777 + 1)3777 + (777)3(777 + 1) 



(n) 4 



(777)2(777 + 1)2 



(77)4 



EV^ / (n*) 1 (r,t) 

V%l,0),(2,l);(l,0),(3,l)],n(6,a) + %l,l),(2,0);(l,l),(3,0)],n(i>,a) 



a,be{0,l} rjtt 



(r,t) , (r,t) \ 

3[(l,0),(2,l);(l,l),(3,0)],n(6,a) + %,l),(2,0);(l,0),(3,l)],n(M) > / 

_ 4(777 + l)|(2m- l) 2 



(108) 



Starting with the first term, indexed by [(1,0), (2, 1); (1,0), (3, 1)], arguing as for the case p = 4, with 
a + /3 = 3 as in (|105j) we obtain the conditioning event (Xqi, -^02, -X03) = (0, e2 + 1, 1 + e 4 ) and may write 



%i,0),(2,l);(l,0),(3,l)],n(6,a) — X] /a.^.»m,m+i 



(109) 



a,&,e{0,l} 



where 



#a,/J = {(ei, e 2 ) : 1 + l(e 2 + 1=0) + l(e a + 1 = 0) = a} = {(e 4 , e 2 ) : 1 + l(e 2 = 1) + l(e x = 1) = a} 
and, with pfe(ci) as in (|107|l . 



Explicitly, the required sets are given by 

£3,0= {(1,1)} £2,1 ={(0,1), (1,0)}, £1,2 = {(0,0)} and £ . 3 = 0, 
and calculating the weights, using ^3(1) = £3(2) = (m + l) 2 (2m — 1)/(t7) 3 , we obtain 

(777 + l)|(2m - l) 2 2(777 + 1)^(2777 - l) 2 , (777 + 1)^(2777 - l) 2 

W3,o = — , 102,1 = and 101,2 = 



[n)l 



Ml 



(")§ 



3G 



Hence, applying Corollary 14. II in (|109[) yields 

4li),(2,l);(l,0),(3,l)],n(6,a) ( 110 ) 

a,b.e{0.1} 

(m + l)|(2m- l) 2 (I 



(n 



1 2 



— (1 + 3A ni i, nm m+1 + 3A„ : 2,n m>m + i 



+ ^■(1 — A ni i, n?Il _ m + 1 — A„.2,n mi „ l+ i + A n! 3 inmm+1 ) 



Similarly, for the term in (jl08[) subscripted by [(1, 0), (2, 1); (1, 1), (3, 0)] the conditioning event becomes 
(Xoi, X02, -^03) = (1, S2 + 1, ei), leading to the collection of sets 

#3,0 = 0, # 2 ,i = {(1,0)}, E ia -{(0,0), (1,1)}, and ^,3 = {(0,1)}, 

that give rise to the same weights as for the terms in (II 10|) . Hence, this term contributes 

(m + l)|(2m-l) 2 (I 

^2 \8 A n,l,n m , m + i - A n,2,n m , m+1 - <\n,3,n m ,m+i) V iiJ J 

+ ^(1 — ^n,l,n,„ >m+ i — A n! 2,n mm + 1 + An i 3.n m ,„ +1 ) 

+ g(l — 3A„ i i ! n m _ m + 1 + 3A rl! 2,n m m + 1 — A nj 3 iIlm m + 1 )^ , 

and adding to (|110[) results in 

(m + l)I(2m-l)Vl fl , 3A \ \ A 

(m + l)|(2m- l) r 



(n)§ V4 V ^ ' —"•-'-»"»+" ' 2 V ^ "»^.^.™+^ 1 4 



(n 



2 • (112) 

3 



For the term in the sum subscripted by [(1, 1), (2, 0); (1, 1), (3, 0)] the conditioning event is (X i, X 2, ^03) = 
(0, e\, 62), resulting in sets 

£3,0 = {(0,0)}, E 2A = {(0,1), (1,0)}, E li2 = {(1,1)} and £ ,3 = 

and a term, therefore, agreeing with (|110p . Lastly, for the term subscripted by [(1, 1), (2, 0); (1, 0), (3, 1)] the 
conditioning event is (X01, X02, -^03) = (1, ei, 1 + e^), the collection of sets becomes, 

#3,0 = 0, #2,1 = {(0,1)}, #1,2 = {(0,0), (1,1)}, and #0,3 = {(1, Q)}, 

and we thus obtain a term agreeing with (|111[) . Hence these last two terms also contribute the factor (|112[) . 
and summing over the cases (r, t) = (to, to + 1) and (r, t) = (to + 1, to) yields (|108[) . 
Lastly, we show that for r, t € {to, to + 1}, r =/= t, 

fl , [(l,0),(2,l);(l,0),(2,l)],n(b,o) = (l + 2 ^n,l,n m , m+ i + A n! 2,n m>m + 1 ) / n2 ^ " (H^) 
a,he{0,l} V >2 

In this case the conditioning event is simply {Xq\, X02) = (0, 0), and selecting one on and one off switch in 
each of the stages r and t yields the weight W2,o = P2(l) 2 = 4(to + 1)2/(^)2; a Pply m g Corollary 14.11 now 
completes the verification of (|113p . 
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